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Aeronautics and Space Administration, Langley Research Center, Hampton, VA, and The 
Strategic Defense Initiative Organization, Key Technologies Office, Washington, D.C. 
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1.0 INTRODUCTION 


1.1 BACKGROUND 

A computer model for atomic oxygen exposure of Earth satellite surfaces which may 
shadow or reflect on one another is described in reference 1. This model is termed the computer 
model for microenvironment atomic oxygen exposure. Included in the suite of programs 
comprising the model are the microenvironment atomic oxygen exposure program SHADOWV2, 
an input file builder for SHADOWV2 called MDDB, and a program to calculate the average atomic 
oxygen exposure to unshielded surfaces called AVESHAD. Also used with these programs is 
TECPLOT, a commercially available graphics program. 

The computer model for microenvironment solar exposure was developed to predict solar 
exposure to satellite surfaces which may shadow or reflect on one another. Included in the new 
model are the microenvironment solar exposure program SOLSHAD and an upgraded version of 
the microenvironment data deck builder MDDB Version 2.0. As in the atomic oxygen model, 
TECPLOT is used for graphics display. The new model has many features in common with or 
similar to those of the computer model for microenvironment atomic oxygen exposure. The 
application of the new microenvironment solar exposure model is described in this document 

1.2 OBJECTIVE 

This document describes the technical features of the computer model for 
microenvironment solar exposure and provides directions for its operation. 

1.3 SCOPE 

Precise treatment of solar reflections requires data on solar spectral reflectance and 
absorption properties of materials used in satellite design. The current model does not explicitly 
include the spectral nature of solar reflectance and absorption. Instead, it assumes these properties 
are constant for a given event and, based on them, calculates solar exposure to surfaces in terms of 
cumulative equivalent Sun hours (CESH); that is, the number of hours of normal incidence solar 
exposure to a hypothetical surface which equals die solar exposure on the actual surface. The 
model assumes a circular satellite orbit around a spherical Earth. This approximation is quite 
accurate for near-circular orbits, but becomes increasingly less accurate as the orbits become more 
highly elliptical. 

A predecessor of SOLSHAD, which is functionally equivalent to the program described in 
this document, was used to model solar exposure to the Long Duration Exposure Facility (LDEF) 
experiment trays (refs. 2 and 3). 
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2.0 PROGRAM OPERATION 


2.1 DESCRIPTION OF COMPUTER MODEL FOR MICROENVIRONMENT 

SOLAR EXPOSURE 

The computer model for microenvironment solar exposure comprises a suite of programs. 
The program which calculates solar exposure and the program to create input files were developed 
at Boeing. The results are displayed using the commercial software package TECPLOT. 
(TECPLOT is a product of Amtec Engineering, Inc., P.O. Box 3633, Bellevue, WA 98009-3633.) 
Elements of the model are shown in the conceptual block diagram of figure 2.1-1; this section 
gives a brief description of these elements. A more detailed description of the model is given in 
section 3.0. 

After the model description, detailed instructions on how to install and operate the various 
pieces of the model are given. They are presented in the order that they will be used: first, how to 
construct an input file for SOLS HAD, including the operation of MDDB; second, how to run 
SOLSHAD; and last, how to view the results. Appendix A contains sample files generated by 
these programs. Appendix B contains the formats of files for TECPLOT. 

2.1.1 Description of Program Modules 

Microenvironment Computer Model for Solar Exposure. SOLSHAD is a 
FORTRAN program which calculates the solar exposure in CESH to satellites in circular orbit 
around the Earth. The satellite is assumed to have a fixed orientation with respect to its motion in 
orbit; that is, each surface maintains a constant orientation with respect to the direction of satellite 
motion and with respect to an axis radially outward from the center of the Earth through the 
satellite. The LDEF satellite is an example of a satellite with this orientation. The total solar 
exposure calculated includes direct solar exposure, solar exposure due to specular and diffuse 
reflection of direct solar exposure from satellite surfaces. Earth-reflected solar exposure, and solar 
exposure due to specular and diffuse reflection of Earth-reflected solar exposure from satellite 
surfaces. During all calculations the satellite is assumed to be in a circular orbit around a spherical 
Earth. Solar exposure for missions or events of arbitrary lengths from fractions of an orbital 
period to many years may be modeled. 

’Die satellite surface is defined by the following geometrical surface elements: trapezoids 
and sections of cylinders, cones, disks, and spheres. Using combinations of these elements, 
complex satellite surfaces may be defined. 

SOLSHAD calculates the solar exposure on satellite surfaces as follows. The satellite 
geometry description is read in and a grid of points is constructed representing each of the surfaces 
on the satellite structure. A visibility matrix is built up by ray tracing from every point on each 
surface to every other point to determine which surfaces may reflect light to themselves or to other 
surfaces. The visibility matrix is used to make ray tracing more efficient later in the program. 

Solar exposure is calculated using a large enough number of Sun positions to give a 
statistically valid representation of all possible Sun positions during the event time. Each Sun 
position is generated by randomly selecting a day of year and time of day and determining the 
subsolar Earth latitude and longitude. Along with each Sun position, a random satellite position is 
chosen from those allowed by the specification of the satellite orbit. The Sun and satellite position 
determine the orientation of die satellite surfaces with respect to die Sun. 
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Surface 

properties 


Figure 2.1-1. SOLSHAD Flow Diagram. 

For each Sun position a ray is traced from each point on the satellite structure to the Sun. 
Those rays which arc unblocked either by other satellite surfaces or by the Earth result in direct 
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solar exposure. In this manner, the direct solar exposure is calculated as the fraction of all possible 
rays from the Sun which reach each point of the structure. Multiplying this fraction by the event 
time in hours gives the CESH of direct solar exposure for each point 

Each ray from the Sun which strikes a point on the structure is Monte Carlo scattered from 
random positions near the point based on the surface properties: specular reflectivity, diffuse 
reflectivity, and absorb tivity. If a scattered ray strikes another surface, it is Monte Carlo scattered 
again. If the scattered ray is absorbed or scattered to space, scattering is terminated for that ray. In 
this manner, CESH are built up for specular and diffuse exposure from direct solar exposure. 

Earth-reflected solar exposure to the satellite structure is handled in a similar way. The 
portion of the Earth which can be seen from the satellite is divided into a grid of points. Each point 
on Earth is tested to see if it is in sunlight. If so, a ray from the Sun to that point is 
deterministically diffusely reflected to die satellite with a diffuse reflectivity of the Earth defined to 
be equal to the average Earth albedo under the Satellite orbit. This Earth-reflected ray is treated just 
as direct solar rays are for exposure on the satellite structure and direct Earth-reflected CESH and 
scattered Earth-reflected CESH are accumulated. 

When all Sun and satellite positions have been considered, the final output of SOLSHAD is 
written to two files: the first file, written to standard output, contains a complete listing of all the 
information generated by SOLSHAD including the CESH from direct and scattered solar and 
Earth-reflected radiation; the second file consists of information that is used by TECPLOT to 
generate 3D graphs of the structure, which display the solar exposure to surfaces as a color 
spectrum. Once displayed in TECPLOT, the structure can be rotated, translated, scaled, and 
printed. One can also draw a line across the surface, and convert CESH along that line into a 2D 
graph of CESH versus distance. 

Microenvironment Data Deck Builder Version 2.0. MDDB Version 2.0 is an 
enhanced and expanded version of the FORTRAN program MDDB Version 1 (ref. 1). While 
MDDB Version 1 allows users to build or edit input files for the microenvironment atomic oxygen 
program SHADOW Version 1.X (ref. 1), MDDB Version 2 allows users to build or edit 
interactively input files for SHADOW Version 1.X, SHADOW Version 2.0 (an unreleased, 
optimized version of SHADOW developed by Boeing), and SOLSHAD, as well as to create any of 
these input files from any of the others. This interconversion is allowed because the three input 
files have large parts in common. The geometry definition sections of all three are identical except 
for different surface properties for the atomic oxygen programs and SOLSHAD. Input for 
SHADOW Version 1.X and Version 2.0 is identical except that the Version 2.0 input file does not 
require a visibility matrix. Calculation control parameters differ between SHADOW and 
SOLSHAD. 

In the remainder of this document, whenever MDDB is mentioned without a version 
number, the reference is to MDDB Version 2.0. 

During the geometry definition phase, MDDB prompts the user for inputs defining surfaces 
or allows the user to modify or delete existing surfaces, to add new surfaces, or to review or 
display the surfaces in a structure as far as it has been built Surface display by the TECPLOT 
plotting package may be selected throughout the geometry definition phase of input file 
construction. All surfaces of a structure are displayed with an arrow pointing outward from the 
active side of each surface. After a surface geometry has been defined, the surface may be 
translated or rotated to reorient it. This allows the user to define the surface in a convenient set of 
coordinates and then move it into its true position in a structure. 
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Because the ray trace routine in SHADOW and SOLSHAD must check every surface to 
determine which if any surface a ray strikes, the ray trace can be sped up by introducing a matrix, 
called the visibility matrix, which specifies which surfaces can see a given surface. Only these 
surfaces need be considered during ray tracing. SHADOW Version 2.0 and SOLSHAD generate 
the visibility matrix automatically; SHADOW Version 1.X requires die user to input the visibility 
matrix. When building SHADOW Version IX input files, MDDB prompts the user to input the 
visibility matrix. The user may review or modify die visibility matrix or display the structure with 
TECPLOT at any time during entry of the visibility matrix. 

It should be noted that at all times during display of the structure, the user is out of MDDB 
(even though MDDB is still opened) and in TECPLOT. Hence, all features of TECPLOT are 
operative: the plot may be rotated, redisplayed as a wireframe or with shading, probed to determine 
the coordinates of a surface, printed, etc. When done viewing die TECPLOT display, the user 
switches back to MDDB while leaving TECPLOT open. 

After completion of geometry definition, MDDB prompts the user for parameters 
controlling calculation and writes the newly created input file to a name of the user's choosing. 


2.1.2 Mainframe and Workstation Platforms 

The microenvironments computer model for solar exposure SOLSHAD and its companion 
program MDDB generally operate on two platforms, a mainframe and a PC, respectively. MDDB 
is designed to operate in a Microsoft Windows Version 3.1 environment in conjunction with the 
commercial software package TECPLOT Version 5.0x. If TECPLOT is available on a more 
powerful workstation (e.g., Sun, Silicon Graphics, Apollo, VAX), the FORTRAN code for 
MDDB can be compiled and run on any of those platforms as well, although some modifications 
may have to be made to the MDDB-generated binary TECPLOT file to be compatible with the 
installed version of TECPLOT. SOLSHAD was developed on a CONVEX C2 supercomputer. It 
generates two files: one that contains all the information from the run, and another that can be 
imported directly into TECPLOT. The TECPLOT file is transferred to the PC or workstation, 
where 2D and 3D plots are generated. 

2.1.3 Input Files 

SOLSHAD Input File. SOLSHAD requires a single input file, which it reads from 
standard input This file contains a parametric description of the satellite geometry and surface 
properties and calculation control parameters. A detailed description of the content of this file is 
given in section 2.2. Most users will find it convenient to use MDDB to construct or edit 
SOLSHAD input files. 

MDDB Version 2.0 Input File. All user input to MDDB is entered interactively from 
the keyboard in response to prompts displayed on the screen. MDDB prompts the user to supply 
input for building new SHADOW or SOLSHAD input files or for editing existing SHADOW or 
SOLSHAD input files. (See sec. 2.2 for a detailed description of the SOLSHAD input file.) 

2.1.4 Output Files 

SOLSHAD Output Files. SOLSHAD produces two output files. The first is directed 
to standard output. This file contains a copy of the SOLSHAD input file, a summary of the object 
geometry, statistics about the ray tracing calculations, a table of CESH at each node, and a listing 
of the connectivity matrix. 
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The second file, named TAPE7 as a default, is an ASCII file that contains the coordinates 
of all the nodes, the CESH at each node, and a connectivity matrix that describes how the nodes 
arc connected to each other. This file is formatted for processing by PREPLOT (a companion 
program to TECPLOT) to create a binary input file for TECPLOT. TECPLOT displays this file as 
a plot of the object with the CESH color-coded on the surfaces. An exact description of this file is 
given in Appendix B. 

MDDB Version 2.0 Output Files. As described in section 2.1.3, MDDB displays 
prompts on the screen for keyboard input The output file produced by MDDB is a SHADOW or 
SOLSHAD input file. MDDB also produces a binary TECPLOT file and an ASCII TECPLOT file 
whenever a TECPLOT display of the object is selected. The binary file, MDDB.PLT, may be used 
directly by the PC version of TECPLOT 5.0x. The ASCII file, MDDB.PRE, must be processed 
by the TECPLOT companion program PREPLOT before a TECPLOT display can be generated. 
The ASCII file is provided so that TECPLOT displays can be made using versions of TECPLOT 
which are incompatible with the binary file. An exact description of these two files is given in 
Appendix B. 

2.2 INPUT FILE CONSTRUCTION 
2.2.1 SOLSHAD Input File 

The SOLSHAD input file divides itself into two natural parts. The first part is a geometry 
description which defines the surfaces comprising the object or structure whose solar exposure is 
to be modeled. The second part is a set of parameters which controls calculations in SOLSHAD. 
The two parts are arranged in consecutive order to construct the input file. Figure 2.2. 1- 1 shows 
an example of a SOLSHAD input file, and figure 2.2. 1-2 shows a symbolic representation of the 
format for the five surfaces. (Experienced users who need only look up the formats of the 
geometry inputs will find this second figure useful.) The SOLSHAD input file must be directed to 
standard input when running SOLSHAD. 

Geometry Description. A structure is constructed from a series of surfaces. These 
surfaces may be any combination of up to 100 trapezoids or segments of cylinders, cones, disks, 
or spheres. It must be emphasized that because they are surfaces and not solids, trapezoids and 
disks have no depth, cylinders and cones are open on the ends, and spheres are hollow. Each 
surface has a front or positive side and a back or negative side. Only one side of a surface may be 
selected for exposure to solar radiation. This side is called the active side. Any solar radiation 
striking the unselected side of a surface will be ignored. This may lead to erroneous results. 
Therefore, the user must be careful when defining which surfaces will be exposed. If the user 
wants to model solar radiation exposure to both sides of a physical surface, he or she must model 
the physical surface with two closely spaced surfaces whose active faces point away from each 
other. 


Each surface in a structure is divided into a grid of nodes. Each surface must have at least 
one node. The total number of nodes for all surfaces may not exceed 5,000. The number of nodes 

on a surface is the product of the number of nodes NC in the xi (£) direction and NN in the eta (rj) 
direction (see record A2 in fig. 2.2. 1-1). Given the choice of constructing a structure having a 
given total number of nodes with many surfaces and a few nodes on each or a few surfaces with 
many nodes on each, the latter choice is to be preferred because the time to trace a ray in 
SOLSHAD is proportional to the number of surfaces. Therefore, the program will execute more 
quickly with only a few surfaces. 
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Record 12345 678 9012345 67 8 901234 5 6789012345 678901234 56789012345678901234 567 8 901234 567 890 


A1 

SAMPLE GEOMETRY WITH 

5 SURFACES 

SOLSHAD 


A2 

15 5 

T T T 




A3 


-4.0000 


-4.0000 15.0000 


A4 


-4.0000 


4.0000 15.0000 


A5 


4.0000 


4.0000 15.0000 


A € 




1.0000 


A7 

PLANE 





A8 

PLANE MATERIAL 




A 9 

0.5000 

0.5000 

0.0000 

SPEC, DIFFUSE, ABSORPTION 


A2 

2 5 24 

T T T 




A3 


0.0000 


0.0000 1.0000 


A4 


0.0000 


0.0000 0.0000 


A5 


1.0000 


0.0000 0.0000 


A 6 


3.0000 


1.6667 

360.0000 

A7 

CYLINDER 





A8 

CYLINDER MATERIAL 




A 9 

0.2500 

0.7000 

0.0500 

SPEC, DIFFUSE, ABSORPTION 


A2 

3 6 24 

T T T 




A3 


0.0000 


0.0000 1.0000 


A4 


0.0000 


0.0000 5.0000 


A5 


1.0000 


0.0000 0.0000 


A € 


3.0000 


1.6667 0.0000 

360.0000 

A7 

CONE 





A8 

CONE MATERIAL 




A 9 

0.3000 

0.6000 

0.1000 

SPEC, DIFFUSE, ABSORPTION 


A2 

4 4 24 

T T T 




A3 


0.0000 


0.0000 1.0000 


A4 


0.0000 


0.0000 0.0000 


A5 


1.0000 


0.0000 0.0000 


A 6 


6.0000 


0.5000 

360.0000 

A7 

DISK 





A8 

DISK MATERIAL 




A9 

0.0500 

0.9000 

0.0500 

SPEC, DIFFUSE, ABSORPTION 


A2 

5 6 24 

T T T 




A3 


0.0000 


0.0000 1.0000 


A4 


0.0000 


0.0000 12.5000 


A5 


1.0000 


0.0000 0.0000 


A € 


2.5000 


0.0000 180.0000 

360.0000 

A7 

SPHERE 





A 8 

SPHERE MATERIAL 




A 9 

0.0500 

0.9000 

0.0500 

SPEC, DIFFUSE, ABSORPTION 


A10 

0 0 0 

F F F 

END OF GEOMETRY 


B1 

1994 6 

19 0 

.0000 

Y M D H START 


B 2 

1994 10 

27 12 

.0000 

Y M D H END 


B 3 

0.0000 

360.0000 


ORBPOS ( 1 ) , ORBPOS { 2 ) 


B 4 

0.0000 

360.0000 

F 

PRECES (1) , PRECES (2) , TIMCOR 


B 5 

28.5000 



ORBINC 


B 6 

0.2460 6. 

3780E+03 

4 . 5000E+02 ALBEDO, RE, HSAT 


B7 

10 1 

10 

1000 TFT 

MAXRAY , MAXALB , NALB, ISUNPOS, SHORTL, CHECK 

, INACT 

B 8 

5 5 



NTHETA, NPHI 


B 9 

997531 


ISEED 


Record 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 


Figure 2.2. 1 - 1 . Sample SOLSHAD Input File. 
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Record 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

Al 

COMMENT ABOUT THE GEOMETRY 




A2 

1 NC NN T T T 

(T T T - DEBUG FRONT 

NODESEL) 


A3 

Plx 

Ply 

Plz 


A4 

P2x 

P2y 

P2z 


AS 

P3x 

P3y 

P3z 


At 



X 


A7 

COMMENT ABOUT THE TRAPEZOID 




A8 

MATERIAL TYPE 




A9 

SPECULAR REFLECTIVITY 

DIFFUSE REFLECTIVITY 

ABSORPTION 


A2 

2 NC NN T T T 

(T T T * DEBUG FRONT 

NODESEL) 


A3 

EX 

Ey 

Ez 


A4 

Ax 

Ay 

Az 


AS 

Pref x 

Pref y 

Pref z 


At 

R 

H/R 


e 

A7 

COMMENT ABOUT THE CYLINDER 




A8 

MATERIAL TYPE 




A9 

SPECULAR REFLECTIVITY 

DIFFUSE REFLECTIVITY 

ABSORPTION 


A2 

3 NC NN T T T 

(T T T = DEBUG FRONT 

NODESEL) 


A3 

Ex 

Ey 

Ez 


A4 

Ax 

Ay 

Az 


AS 

Pref x 

Pref y 

Pref z 


At 

Rmax 

H/Rmax Rmin/Rmax 

e 

Al 

COMMENT ABOUT THE CONE 




A8 

MATERIAL TYPE 




A9 

SPECULAR REFLECTIVITY 

DIFFUSE REFLECTIVITY 

ABSORPTION 


A2 

4 NC NN T T T 

(T T T = DEBUG FRONT 

NODESEL) 


A3 

Ex 

Ey 

Ez 


A4 

Ax 

Ay 

Az 


A 5 

Pref x 

Pref y 

Pref z 


At 

Rout 


Rin/Rout 

e 

Al 

COMMENT ABOUT THE DISK 




A8 

MATERIAL TYPE 




A9 

SPECULAR REFLECTIVITY 

DIFFUSE REFLECTIVITY 

ABSORPTION 


A2 

5 NC NN T T T 

(T T T = DEBUG FRONT 

NODESEL) 


A3 

Ex 

Ey 

Ez 


A4 

Ax 

Ay 

Az 


AS 

Pref x 

Pref y 

Pref z 


At 

R 

♦1 

$2 

e 

Al 

COMMENT ABOUT THE SPHERE 




A8 

MATERIAL TYPE 




A9 

SPECULAR REFLECTIVITY 

DIFFUSE REFLECTIVITY 

ABSORPTION 


A10 

0 0 0 F F F END OF 

GEOMETRY 




B 1 
B 2 
B 3 
B 4 
B 5 
B 6 
B 7 
B 8 
B 9 

Record 


IYR{1) IMO(l) IDM(l) GMT (1 ) 
IYR (2 ) IMO ( 2 ) IDM ( 2 ) GMT ( 2 ) 
ORBPOS (1 ) ORBPOS ( 2 ) 

PRECES(l) , PRECES{2) ,TIMCOR 
ORBINC 

RE HSAT 

MAXALB NALB ISUNPOS 
NPHI 


OF EVENT START 
OF EVENT END 

RANGE OF SATELLITE POSITIONS ON ORBIT 

NODAL CROSS LONG, PRECESSION RANGE, TIME CORR (T,F) 


SHORTL (T) CHECK (F) INACT 


ALBEDO 
MAXRAY 
NTHETA 
ISEED 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 
NOTE: All entries are free format except for record A2 which is format (I2,2I5,3L2) and records A3 through A6 
which are format (4F20.10), Some entries in these records are blank. 


Figure 2.2. 1-2. Symbolic SOLSHAD Input File. 
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Solar exposure is calculated in each node on the structure. The resolution of detail of 
exposure on a surface thus depends on the number of nodes on a surface. The user should inspect 
his structure carefully to select noding which gives high resolution in areas where exposure may be 
expected to vary greatly and lower resolution where the exposure is expected to be more slowly 
varying. 

Notwithstanding the foregoing, the grid of nodes should always be fine enough so that the 
grid points approximate the shape of the surface reasonably well. This is because the routine 
which derives the visibility matrix does so by ray tracing from each grid point to every other grid 
point in the structure. If the surfaces are not well approximated by the grid, an incorrect visibility 
matrix may result. The following suggestions for minimum gridding of surfaces may be helpful. 
Trapezoids and small portions of other surfaces which are approximately plane quadrilaterals 
should have at least one node which defines the four comers of the quadrilateral. The grid of 
nodes for sections of disks should be fine enough to approximate die perimeter of the disk. For 
nonplanar surfaces such as sections of cylinders, cones, and spheres, the grid should be fine 
enough to give a reasonable representation of the three-dimensional surface. 

Because surfaces in SOLS HAD are not solids, they may overlap or intersect in ways that 
solids cannot In such cases, a surface may be exposed to solar radiation on one side of the 
intersection or overlap and shielded from it on the other. In some cases the user may wish to take 
advantage of this; in others, surfaces may inadvertently intersect or overlap, causing unanticipated 
results. An example of the latter is when two trapezoids are intended to have a common edge; but 
due to numerical inaccuracy or roundoff, the two trapezoids intersect or overlap near their intended 
common edge. Whether intended or not, points on the nodes of the intersecting or overlapping 
surfaces will often show anomalous (and incorrect) solar exposure. Unintended intersections and 
overlaps are especially easy to generate when triangles are approximated by trapezoids with a very 
small (but not zero) ratio of the short-to-long parallel side to avoid trouble with TECPLOT (see 
Possible Problems at the end of sec. 2.2.3) When unintended intersections or overlaps occur, the 
user should adjust the geometry of surfaces in the structure to remove the problem. It should be 
noted that small gaps between the edges of surfaces seldom have a great effect on modeled solar 
exposures. When the intersections or overlaps are intended, the region of incorrect solar exposure 
may be minimized by using a fine grid. 

Coordinate Axes Orientation. Coordinates defining surfaces in a structure are defined in a 
single, well specified, orthogonal, right-hand Cartesian coordinate system. In this coordinate 
system (see fig. 2.2. 1-3), the origin is at the satellite. The positive x axis is pointed radially 
outward along a line from the center of the Earth to space. The positive z axis is pointed 
perpendicular to the x axis in the direction of satellite motion along its orbit. The positive y axis is 
oriented so the (x,y,z) form a right-handed coordinate system. Users building structures with 
MDDB will often find it easiest to build structures in a convenient user-defined coordinate system 
and then rotate the entire structure to the coordinate system required by SOLS HAD. 

Refer to figures 2.2. 1-1 and 2.2. 1-2 for samples of SOLSHAD input files. The geometry 
description takes place in records labeled A1 through A10 for this example with five surfaces. 
Before giving complete definitions for these records, a few words about the structure of the surface 
definition records are in order. Record A1 is a text description of the structure and is entered only 
once. Records A2 through A9 are repeated once for each surface in the structure. Records A2 
through A6 define the surface type and geometry. Record A7 contains a text description of the 
surface. Records A8 and A9 describe and define surface properties. Record A10 terminates the 
geometry input. 
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Z (HEADING 
DIRECTION) 


Y 


X (SPACE) 



Figure 2.2. 1-3. Geometry Coordinate Definition. 


Record A1 Format (A) 

Entry 

Description of geometry. 

This record is entered once only. 


Columns 

1-60 


T ype Variable 

character" HEADER 


10 





Records A2 through A9 are repeated in sequence for each surface. 


1-2 


T 


integer 


Variable 


IP(IPS,1) 



Record A2 Format (12,215, 3L2) . 1 1 

— I Columns 1 ~ 


Surface type 

1 = trapezoid 

2 = cylinder 

3 = cone 

4 = disk 

5 = sphere. 

Number of nodes in xi direction. Refer to figures 
2.2. 1-4 through 2.2. 1-8 for orientation of nodes. 


Number of nodes in eta direction. Refer to figures 
2.2. 1-4 through 2.2. 1-8 for orientation of nodes. 


Debug flag. If T (true), select normal output 1 . If LASTF 

F (false), select full debug output for geometry. F 
(false) is rarely used because it generates an 

excessive ly large output . — . — . . _>vv t 

Flag to designate active side of surface. If T (true), 15-16 logical F(15)ot 

select front (positive) side of surface. If F (false), rKUlN 1 

select back (negative) side of surface . . — »i 4 

Flag to select whether CESH on this surface is to 17-18 logical b 6) or 

be calculated. If T (true), calculate CESH on NODESEL 

surface 2 . If F (false), do not calculate CESH on 
surface. F (false) is rarely used because careless 

use may cause erroneous results. | 1 1 

The contents of records A3 through A6 depend on which surface has been selected in 
record A2. 

Surface Type 1 • TRAPEZOID. The trapezoid (fig. 2.2. 1-4) is defined by three points and 

the ratio of the lengths of the parallel sides. This ratio, X, must be 0 ^ X ^ 1. Let PI and P2 be the 
points located at the ends of the longer of the two parallel sides. P3 is an end point on the shorter 

parallel side such that P2 is connected to PI and P3. A triangle can be made by setting X = 0. The 
positive side of the surface is defined such that the comer points are ordered in counterclockwise 
manner when the positive side is visible. This is the conventional right-hand rule. 


Record A3 Format (3F20.10) 




Variable 



'MDDB automati call y selects this entry. Users wishing to enter F (false) must edit their input files outside MDDB. 
2 MDDB automatically selects this entry. Users wishing to enter F (false) must edit their input files outside MDDB. 
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Record A4 Format (3F20.10) 


Entry 

HS!!IT!nTnH 

Type 


X coordinate of P2. 

1-20 

real 

■a imi 

Y coordinate of P2. 

21-40 

real 


Z coordinate of P2. 


real 

wmmsm 


Record AS Format (3F20.10) 


Entry 


Type 

THlf 

X coordinate of P3. H 

1-20 

real 


Y coordinate of P$. 


real 


Z coordinate of P3. 

4TroO 

real 

P(IPS,9) 


Record A6 Format (20X,F20.10) 


Entry 

msEsmsi 

Type 

■TOM 

Ratio of short to long parallel sides of trapezoid, X. 

21-40 

real 

mssm 


P4 P3 



Figure 2.2. 1-4. Trapezoid. 


Surface Type 2: CYLINDER. The cylinder (fig. 2.2. 1-5) is defined by an axial vector E, a 
center point A, a reference point P re f, the radius R, the ratio of height to the radius H/R, and the 

angular extent 0. The axial vector E points from the center point A through the cylinder. The 
reference point P re f and the axial vector E specify a half plane from which the angular extent of the 

cylinder 0 is measured. 0 is measured in the right hand-rule sense with respect to the direction of 
E. The positive side of the cylinder is the outside. 
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Record A3 Format (3F20.10) 


Entry 

Columns 

Type 

Variable 

X coordinate of axial vector E. 


real 


Y coordinate of axial vector E. 

mriZL'wm 

real 

P(IPS,2) 

Z coordinate of axial vector E. 

41-60 

real 

P(IPS,3) 


Record A4 Format (3F20.10) 


Entry 

Columns 

Type 

Variable 

X coordinate of center A. 


real 

P(IPS,4) 

Y coordinate of center A. 

— in i 

real 

msmsam 

Z coordinate of center A. 

WMSBSm 

real 

P(IPS,6) | 


Record A5 Format (3F20.10) 


Entry 

Columns 

Type 

Variable 

X coordinate of reference point P re f- 


real 


Y coordinate of reference point P rc f- 

21-40 

real 

P(IPS,o) | 

Z coordinate of reference point Pref- 

41-60 

real 

■aaaaaM 


Record A6 Format (2F20.10,20X,F20.10) 


Entry 

Columns 

Type 

Variable 

Radius of cylinder R. 

1-20 

real 

C(2) 

Ratio of cylinder height to radius H/R. I 

21-40 

real 

C(3) 

Angular extent of cylinder 0. 0 < 0 Ss 360 degrees. 

61-80 

real 

C(5) 


Surface Tvne 3: CONE. The cone (fig. 2.2. 1-6) is defined similarly to the cylinder. The 
cone is defined by an axial vector E, a center point A, a reference point Pref, a maximum radius 
Rrnax, the ratio of height to the maximum radius H/Rmax, the ratio of the minimum to maximum 

radius Rmin/Rmax, and the angular extent 0. The axial vector E points from the center point A 
through the cone. The reference point P re f and the axial vector E specify a half plane from which 
the angular extent of the cone 0 is measured. 0 is measured in die right-hand rule sense with 
respect to the direction of E. The positive side of the cone is die outside. 


Record A3 Format (3F20.10) 


Entry 


Type 

Variable 

X coordinate of axial vector E. 


real 

P(IPS,1) 

Y coordinate of axial vector E. 

21-40 

real 

P(IPS,2j 

Z coordinate of axial vector E. 

41-60 

real 



Record A4 Format (3F20.10) 


Entry 

Columns 

Type 

Variable 

X coordinate of center A. 

r~ 1-20 

real 


Y coordinate of center A. 

21-40 

real 


Z coordinate of center A. 

41-60 

real 
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Figure 2.2.1-5. Cylinder. 


Record A5 Format (3F20.10) 


Entry 

Columns 

Type 


X coordinate of reference point P re f. 

1-20 

real 

■aiasMi 

Y coordinate of reference point P re f. 

21-40 

real 

■a — 

Z coordinate of reference point P re f. 


real 

mssmm 


Record A6 Format (4F20.10) 


Entry j 

Columns 

Type 

Variable 

Maximum radius of cone Rmax- 

T i 'tT * 

real 

wms®m 

Ratio of cone height to radius H/Rmax- 


real 




real 

mm 

Angular extent of cone 0. 0 < 0 < 360 degrees. 

61-80 

real 

KM 
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Surface Type 4: DISK. The disk (fig. 2.2. 1-7) is defined by an axial vector E, a center 
point A, a reference point P re f, an outer radius Rout> 311 inner radius Rjn, the ratio of the inner to 

outer radius Rin/Rout. and the angular extent 0. The reference point P re f and the axial vector E 

specify a half plane from which the angular extent of the disk 0 is measured. Note that 0 is 
measured in the left-hand rule sense with respect to the direction of E. The positive side of the 
disk is the side whose outward surface normal is parallel to E. 


Record A3 Format (3F20.10) 


Entry 

Columns 

Type 

Variable 

X coordinate of axial vector E. 

1-20 

real 

P(IPS,1) 

Y coordinate of axial vector E. 

3130 

real 

P(IPS,2) 

Z coordinate of axial vector E. 


real 

P(IPS,3) 


Record A4 Format (3F20.10) 


Entry 


Type 

Variable 

X coordinate of center A. 

1-20 

real 

P(IPS,4) 

Y coordinate of center A. 

21-40 

real 

P(IPS,5) 

Z coordinate of center A. 

41-66 " 

real 

P(IPS,6) 


Record A5 Format (3F20.10) 


Entry 

Columns 

Type 

Variable | 

X coordinate of reference point P re f. 

1-20 

real 


Y coordinate of reference point P re f- 

21-40 

real 

P(IPS,8) 

Z coordinate of reference point P re f. 

41-60 

real 

P(IPS,9) 
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Record A6 Format (F20.10,20X,2F20.10) 



Columns 

Type 


Outside radius of disk Ro Ut . 

1-20 

real 

HB 1 

Ratio of inside to outside radius of disk Rin/Rout- 

HK&sIH 

real 

C(4) 

Angular extent of cone 0. 0 < 0 < 360 degrees. 
Note that this is measured in the left-hand rule 
sense. 

61-80 

real 

■fl 



Positive 

side 


Surface Type 5: SPHERE. The sphere (fig. 2.2. 1-8) is defined by an axial vector E along 
the polar axis, a center point A, a reference point P re f, the radius R, the azimuthal angular extent 0, 

and two angles defining the polar (zenith) angular extent <)>i and <I>2, with <J>i being the smaller 
angle. The axial vector E points from the center point A through the north pole of the sphere. The 
reference point P re f and the axial vector E specify a half plane from which the azimuthal angular 

extent of the sphere 0 is measured. 0 is measured in the right-hand rule sense with respect to the 
direction of E. The positive side of the sphere is the outside. 
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Figure 2.2. 1-8. Sphere. 


Record A3 Format (3F20.10) 


Entry 

Columns 

Type 

Variable 

X coordinate of axial vector E. 

1-20 

real 

C(l) 

Y coordinate of axial vector E. 

21-40 

real 


Z coordinate of axial vector E. 

■KU&H 

real 

wnsswm 
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Record A4 Format (3F20.10) 


Entry 

■K.HUnl.l-M 

Type 

— wasgtflqMi 

X coordinate of center A. 

P20 

real 


Y coordinate of center A. 

2T35 

real 

C(6) 

Z coordinate of center A. 

41-60 

real 

mtssimm 


Record AS Format (3F20.10) 


Entry 


_Type 


X coordinate of reference point P re f. 

WKME&MM 

real 

KH 

Y coordinate of reference point P re f. 

21-40 

real 

msmm 

Z coordinate of reference point P re f. 

41-60 

real 

mmssm 


Record A6 Format (4F20.10) 


Entry 


Type 

Variable | 

Radius of sphere R. 


real 

■ HIM 

Smaller polar angle limit 4»i. 0 < <(>i < <j >2 degrees. 

| 

real 

C(14) 

Larger polar angle limit <j> 2 . <j>i < <j >2 ^ 180 degrees. 


real 


Azimuthal extent of cone 0. 0 < 0 < 360 degrees. 

61-80 

real 



Record A7 Format (A) 


Entry 

Columns 

Type 

Variable 

A comment about the surface. 

m 

character 

C(31)-C(40) 

or 

P(IPS,17)- 

P(IPS,26) 


Record A8 Format (A) 


Entry 

HSRlRKfTTI^H 


The name of the surface material. 

1-80 1 



Record A9 Free Format 


Entry 

Type 

Variable 

Specular reflectivity of surface for solar radiation. 

real 

SPROP(l,IPS) 

Diffuse reflectivity of surface for solar radiation. 

real 

SPROP(2,IPS) 

Absorptivity solar radiation. 

real 

SPROP(3,IPS) 


The surface properties specify what fraction of incident solar radiation on the surface is 
specularly reflected, diffusely reflected, and absorbed. These three fractions must sum to one. 
Table 2.2.1-1 contains surface properties for Fluorinated Ethylene Propylene (FEP) and aluminum 
used for modeling several cases. 
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Table 2.2.1-1 

Solar Surface Properties of Some Materials 


Material 




Used to Model 

FEP 

0.83 

0.1 

0.07 

2, 3, 4, 5 

Aluminum 



1 

1. 2, 3. 4, 5 


Used to Model: 

1. CRES 

2. LDEF M0001 NRL cosmic ray experiment 

3 LDEF S0069 row 9A thermal control surfaces experiment 

4. LDEF FEP blanket fold tray C5 near row 6 edge 

5 . LDEF FEP blanket fold tray D 1 near row 6 edge 


Entry 

Columns 

SE 

Variable 

0 

1-2 


IP(IPS,1) 

0 

3-7 



"0 

8-12 

integer 

IP(IPS,3) or 
NN 

F (false). 

HUifl 

logical 

F(13) or 
LASTF 

F (false). 

mam 

logical 

F(15) or 
FRONT 

F (false). 

17-18 

logical 



proper geometry termination and for compatibility with MDDB. 

Calculation Control Parameters. The calculation control parameters on the following 
nine, records determine the conditions under which solar exposure takes place. 

The following two records define the beginning and end of the event for which solar 
exposure is to be calculated. 


Record B1 Free Format 


Entry 

Type 

Variable 

Two or four digit year of event start If two digits 
are entered, 0 to 49 is assumed to be 2000 to 2049; 
50 to 99 is assumed to be 1950 to 1999. The 
earliest allowed year is 1950. Four digits are 
always the actual year. 

integer 

IYR(l) 

Month (1 to 12) of event start. 


IMOd) 

Day of month (1 to 31) of event start 

wmsmm 

IDM(l) 

Greenwich Mean Time (GMT) hour and fraction of 
event start. 

real 

GMT(l) 
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Record B2 Free Format 



Type 

Variable 

Two or four digit year of event end. If two digits 
are entered, 0 to 49 is assumed to be 2000 to 2049; 
50 to 99 is assumed to be 1950 to 1999. The 
earliest allowed year is 1950. Four digits are 
always the actual year. 

integer 

IYR(2) 

Month (1 to 12) of event end 


IMO(2) 

Day of month (1 to 31) of event end. 


IDM(2) 

GMT hour and fraction of event end. 

real 

<jMT(2) "1 


The next record specifies the angular range on orbit which should be considered. The 
range is measured from die ascending node crossing. To specify all positions on orbit, set 
ORBINC(l) to 0 degrees and ORBINC(2) to 360 degrees. A smaller angular range on orbit, 
which, for example, might be required for a very short event time, can be specified by specifying 
a restricted range for ORBINC(l) and ORBINC(2). If ORBINC(l) and ORBINC(2) are equal, 
the satellite is fixed at a single position on orbit relative to the ascending node crossing. Through 
the entries on records B3 and B4, the user has great freedom to specify the satellite position on 
orbit during solar exposure. 

Record B3 Free Format 



Type 

Variable 

Start of range of satellite positions measured in 
degrees from ascending node crossing. 

real 

ORBPOS(l) 

End of range of satellite positions measured in 
degrees from ascending node crossing. 

real 

ORBPOS(2) 


The next record specifies the satellite ascending node crossing at the start of the event and 
the range which the orbit may process from there. The use of the variables in record B4 is 
illustrated in the following examples. If it is known that the nodal crossing is at 20 degrees East 
longitude at the start of the mission and processes 30 degrees east to 50 degrees East longitude, the 
user would enter PRECES(l) = 20 and PRECES(2) = 30. Note that PRECES(2) is positive for 
precession toward the east and negative for precession toward die west PRECES(2) may be 
greater than 360 degrees. If PRECES(2) = 0, the orbit does not process. For temporally short 
events the orbital precession is correlated with (proportional to) the time after event start. In such a 
case, TIMCOR is set to T (true). For longer missions (for example, many months with the orbit 
processing about the Earth a number of times), the relationship between ascending node crossing 
longitude and the time after event start is neither obvious nor important In such a case, TIMCOR 
is set to F (false) and the amount of precession is chosen randomly independent of event time. 

This selection gives equal probability of ascending node crossing for any value allowed in the 
range of orbit precession. The standard values of the variables in record B4 for long missions are 
PRECES(l) = 0, PRECES(2) = 360, and TIMCOR = F. This selection gives equal probability of 
the satellite being over any longitude on Earth. 
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Entry 

Type 

Variable 

Initial (start of event) ascending node crossing 
longitude (0 to 360 degrees East longitude). 

real 

PRECES(l) 

Orbital precession (degrees East longitude) allowed 
from initial ascending node crossing). 

real 


T (true) if orbital precession is correlated to event 
time. 

F (false) if orbital precession is not correlated to 
event time. 

logical 

TIMCOR 


Entry 

Type 

Variable 

Orbit inclination (0 to 180 degrees). Values 
between 0 and 90 degrees indicate easterly heading 
orbits and values between 90 and 180 degrees 
indicate westerly heading orbits. 

real 

ORB INC 


Record B6 Free Format 


Entry 

Type 

Variable 

Earth albedo for reflected solar radiation. Earth 
albedo is the average diffuse reflectivity of the 
portion of the Earth under the satellite orbit If 
ALBEDO = 0, no calculation of Earth-reflected 
radiation is performed. Values of albedo range 
from 0 to 1, with 0.246 being typical for 28.5- 
degree orbit inclination. 

real 

ALBEDO 

Average radius of the Earth (km). 6378 km is a 
nominal value. 

im 

RE 

Circular orbit average satellite altitude (km) above 
the surface of the Earth. 

real 

HSAT 


The following two records specify parameters which determine die accuracy of the solar 
exposure calculation and the execution time of SOLSHAD. Sun positions are chosen at random 
times during the event. Assuming a normal distribution of the error in total solar exposure as a 
function of number of Sun positions, the average error in calculated solar exposure to a surface is 
proportional to the square root of the number of Sun positions ISUNPOS. Program execution time 
is proportional to the number of Sun positions. Typically, the user should use 1000 or more Sun 
positions to give a statistically accurate estimate of exposure for long missions. 

No hard and fast rules may be given for specifying the variables on records B7 and B8. 

The user must examine his structure's geometry and use his experience to select values for the 
variables. The following guidelines for long missions may prove helpful. Surfaces which face 
space do not receive Earth-reflected radiation unless it is reflected from another surface. There is 
no need to calculate Earth-reflected radiation for surfaces which cannot receive it. Surfaces which 
face the Earth typically receive several times as much exposure from Earth-reflected solar radiation 
as from direct solar radiation. Vertical surfaces typically receive several times as much direct solar 
radiation as Earth-reflected radiation. The numbers of Monte Carlo scattered rays for direct and for 
Earth-reflected solar radiation should be chosen to reflect the anticipated solar exposure. CESH 
from scattering (specular and diffuse reflection and absorption) is based on constant ray areal 
density; that is, a constant number of rays per unit area. For purposes of scattering calculations. 
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the unit of area is taken to be the average of the nodal areas associated with each grid point in the 
structure. MAXRAY specifies the number of direct solar rays to be scattered per average area and 
MAXALB specifies the same for Earth-reflected rays. The number of rays scattered from the area 
associated with a point equals MAXRAY or MAXALB times the projected area divided by the 
average area. The projected area is the area times the magnitude of die cosine of the angle between 
the ray direction and die surface normal of the area. 

MAXRAY and MAXALB should be set so that several rays are scattered from most nodes 
on the structure. Because the areas associated with nodes on the structure ma y vary by well more 
than an order of magnitude for many typical situations (for example, the areas associated with 
nodes near the vertex of a cone may be much smaller than those near the base of the cone), 
judgment is needed in selecting MAXRAY and MAXALB. A useful way to test whether 
MAXRAY and MAXALB are appropriate is to submit a SOLSHAD job with trial values. Stop the 
job after a minute or two (running the job to completion may take much longer) and examine the 
output for the number of nodes which have no rays scattered from them. If this number is more 
than a few percent of the total number of nodes, larger values of MAXRAY and MAXALB are 
probably in order. Typical starting values of MAXRAY and MAXALB might be 10 and 5, 
respectively. 

The portion of the Earth visible from the satellite is divided into a grid of equal angular 
increments in azimuth angle and in polar angle relative to a line from the satellite to the center of the 
Earth. Record B8 defines this grid. Solar radiation is diffusely reflected from each point on the 
grid to every point on the structure (satellite). Because the time to calculate the Earth-reflected 
radiation exposure from each point on Earth to structure surfaces is approximately as long as the 
time to calculate the direct solar exposure, one wants to minimize the calculation of Earth-reflected 
radiation without introducing significant error into the calculation. Toward this end, the user may 
specify that the Earth-reflected radiation is to be calculated only after a specified number of Sun 
positions. For geometries which receive only a small fraction of their total solar exposure from 
Earth-reflected radiation, the user might set NTHETA and NPHI on record B8 to 5 and NALB on 
record B7 to 10 (calculate Earth-reflected radiation after every tenth Sun position). For an Earth 
facing geometry, the user might set NTHETA and NPHI to 10 and NALB to 1 to calculate Earth- 
reflected radiation at every Sun position. 
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Record B7 Free Format 


Entry 

Number of direct solar rays per average node area 
to scatter during Monte Carlo scattering. This areal 
density of rays is scattered from each node in the 
structure for each Sun position. If MAXRAY = 0, 
only direct solar illumination to the structure is 


integer 


Variable 

MAXRAY 


calculated. 

Solar rays are diffusely reflected from each of the 
NTHETA by NPHI grid points on Earth (record 
B8) to each node on the structure. MAXALB is the 
number of Earth-reflected rays per average node 

area scattered during Monte Carlo scattering. This 

areal density of rays is scattered from each node in 
the structure. If MAXALB =0, only radiation 
reflected from the Earth onto the structure is 
calculated and no Monte Carlo scattering of this 


integer 


MAXALB 


radiation is performed. 

Earth-reflected radiation scattered to the structure 
calculated once every NALB Sun positions. If 
NALB > ISUNPOS, no Earth-reflected radiation 


is 


integer 


NALB 


calculation is done. 

Total number of Sun positions for which solar 

radiation exposure is calculated. 

If T (true), do not list the coordinates of each node 
on the surface of the structure at each step of the 
ray tracing 3 . 

If F (false), list coordinates of each point on the 
surface of the structure at each step of the ray 
tracing. This option should be used only for 
debugging; it produces a very lengthy output. 

If T (true), calculate coordinates of points on ~ 

surfaces of structure, but do not calculate solar 
exposure. 

If F (false), calculate solar exposure on surfaces of 


integer 

logical 


logical 


ISUNPOS 

SHORTL 


CHECK 


structure 4 . 

If T (true), do not print out rays which intersect the 
inactive side of surfaces during the automatic 


logical 


INACT 


generation of the visibility matrix. 
If F (false), print out these rays. 


3 MDDB automatically selects this entry. Users wishing to enter F (false) must edit their input files outside MDDB. 
4 MDDB automatically selects this entry. Users wishing to enter T (true) must edit their input files outside MDDB. 
The true entry has been rendered obsolete by the development of MDDB. 
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Record B8 Free Format 



Type 

Variable 

Number of polar angle increments for the grid on 
the part of the Earth seen by the satellite. 

integer 

NTHETA 

Number of azimuth angle increments for the grid 
on the part of the Earth seen by the satellite. 

integer 

NPHI 


Record B9 Free Format 



Type 

Variable 

Randan number generator seed (an integer 
between 1 and 7 digits). 

integer 

ISEED 


2.2.2 Installation of the Geometry Builder 

The Geometry Builder consists of two separate programs: MDDB (for Microenvironment 
Data Deck Builder) and TECPLOT (a commercial 3D plotting application). These two programs 
arc run simultaneously, using a windowing system that allows switching between them. MDDB 
collects information to build a SOLS HAD or SHADOW input file. At various times it is possible 
to pause in the middle of MDDB, have it write a file for TECPLOT, and then switch to TECPLOT 
in order to view the progress. 

This section describes how to install and run these two programs on a PC in the Microsoft 
Windows environment (version 3.0 or higher). TECPLOT is also available for other types of 
workstations, and MDDB and TECPLOT can be run in a similar manner on these workstations. 
FORTRAN code for MDDB is included; compile it on the workstation and run the two programs in 
separate windows. To recompile the code for Microsoft Windows on a PC, use the Microsoft 
FORTRAN compiler (version 5.0 or higher) with the Microsoft Windows compiler option 
(fl /FeMDDB2 /MW *.for). 

Note : Users of the microenvironment atomic oxygen exposure program SHADOW (ref 1) 
will note that MDDB Version 1 and MDDB Version 2.0 share a number of subroutine names in 
common. Although their names are the same, these subroutines are different and cannot be 
interchanged between the two codes. 

In this section, commands that the user types are written in bold, computer displays are 
written in helvetica, and file names are written in courier. 

Installation for Microsoft Windows 

1. Copy the following files from disk into the same directory that contains TECPLOT: 
mddb2 . exe 
p.mcr 
table .mcr 
shadow. mcr 

Note: The TECPLOT macro files (those ending in . mcr) that you should copy will 
depend on the version of TECPLOT that you are using. For version 5 on a PC, copy the 
files from the directory tec5pc. For version 5 on a workstation, copy the files from the 
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directory tec5ws. For version 6 on a PC, copy the files from the directory tec6pc. For 
version 6 on a workstation, copy the files from the directory tec6ws. 

2. Copy the following file from disk into the same directory that contains Windows: 

tecplot .pif 

3. Edit the TECPLOT configuration file (tecplot . cf g in the directory that contains 

TECPLOT), and make sure that the 3D aspect ratio limit is 100,000 or higher. This will 
prevent flat figures from becoming distorted. 

4. Start up Microsoft Windows by typing win . 

5. Create a Microenvironments Group. 

a. Open up the Program Manager Window. 

b. Choose New from the File Menu. 

c. Choose Program Group and OK. 

d. In the Description box, type Microenvironments. 

e. In the Group File box, type mddb.grp. Hit OK. 

6. Add an icon for MDDB. (Except for the Description box, the values typed in are not case 

sensitive.) 

a. Choose New from the File Menu. 

b. Choose Program Item and OK. 

c. In the Description box, type MDDB. 

d. In the Command Line box, type mddb2.exe. 

e. In the Working Directory box, type C:\tecplot (or whatever the directory is that 
contains TECPLOT and mddb . exe). 

f. Leave “None” in the Shortcut Key box type. 

g. Choose an icon by clicking on the Change Icon button. We suggest a simple 
window icon. 

h. Click on OK. The new icon labeled MDDB should appear in the new window. 

9. Add an icon for TECPLOT. (Except for the Description box, the values typed in are not case 
sensitive.) 

a. Choose New from the File Menu. 

b. Choose Program Item and OK. 

c. In the Description box, type Tecplot. 

d. In the Command box, type c:\windows\tecp!ot.pif (or whatever path gets into 
the windows directory). 

e. In the Working Directory box, type, C:\tecplot (or whatever the directory is that 
contains TECPLOT and mddb . exe). 

f. Leave “None” in the Shortcut Key box. 

g. Choose an icon by pressing the Change Icon button. We suggest a 3D graph icon 
if available. 

h. Click on OK. The new icon labeled Tecplot should appear in the window. Size and 
move the window by clicking and dragging on the bottom comer and top bar, respectively. 

2.2.3 Constructing Input Geometries Using MDDB and TECPLOT 

This section describes how to operate MDDB and TECPLOT simultaneously on Microsoft 
Windows so that an input file can be viewed as it is being constructed. First, it is necessary to 
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open two windows: one that contains MDDB, and the other containing TECPLOT. At various 
times during MDDB operation, you will be asked if you want a TECPLOT display. If you 
respond “yes,” MDDB will generate a file for TECPLOT, and then pause. You will then switch to 
the TECPLOT window, run a macro to call up the file that was generated and display it You may 
then rotate and view the object using standard TECPLOT commands. Once you are finished, you 
will return to the MDDB window, 

MDDB will ask for data with which to build the input file. MDDB version 2 has an 
improvement over version 1 in that it will display default data and ask if you want to keep this data, 
or enter new data. By entering an equal sign (=), the default data will be used; any other input 
will be taken as new data. 

Setting up Windows for TECPLOT and MDDB 

1. Start up Microsoft Windows. 

2. Double-click on Tecplot. After about 10 seconds, the TECPLOT screen will appear. Hit the 
spacebar. 

3. Type <CntlxEscape>. Windows will return, and a box will appear. Double-click on 

Program Manager. 

4. Double-click on MDDB icon. A window will appear labeled MDDB - [Unit *J. 

Constructing Geometries Using MDDB 

1. You are asked for the file format of the file you will be creating. Choose 3 for SOLSHAD. 

(You also have the option of choosing 1 for SHADOW version 1.x input files, or 2 for 
SHADOW version 2.x) 

2. You are asked whether to create a new file or edit an existing one. Enter 1 for a new file, or 2 
for the existing file. If you have modeled a geometry for atomic oxygen exposure using 
SHADOW, and now want to model the solar exposure, you can edit your previous SHADOW 
input file and save it as a SOLSHAD file. 

3. If a new file is chosen, you are asked for an event description, start and end dates, a mission 
file name, and a descriptive header for the geometry. 

4. If an existing file is chosen, enter the file name at the next prompt Then enter the file type: 1 
for SHADOW version 1.x, 2 for SHADOW version 2.x, or 3 for SOLSHAD (if you give the 
wrong file type, the program will not detect the problem, and a program error will occur which 
will halt MDDB). You are then shown the current geometry description and asked to enter a 
new one (or enter = to accept the current one). 

5. You are then presented with the main menu. (If creating a new file, you are automatically sent 
into option 1.) 

SELECT AN ACTION BY NUMBER FROM THE LIST 

1 ADD ANOTHER SURFACE AT END OF SURFACES 

2 MODIFY AN EXISTING SURFACE 

3 DELETE A SURFACE 

4 INSERT A NEW SURFACE 

5 REVIEW SURFACES 

6 COPY SURFACE AND INSERT IT AT A NEW POSITION 

7 EXIT SURFACE OPERATIONS (DONE) 
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6. To add another surface at the end of the existing surfaces: 

<L* Select 1 • 

b. MDDB gives you the “current definition” of the new surface, which is meaningless. 
Press return to enter the new data. 

c. Choose a surface type (1 = trapezoid, 2 = cylinder, 3 = cone, 4=disk, and 5 = sphere.) 

d. You will be prompted for all the information that describes the geometry. (These values 

are described in detail in sec. 2.2.1.) If you are prompted for several values on one 
line, separate them by spaces or commas. For example, to enter the coordinates (4.0, - 
2, 1.3), type 4.0 -2 1.3 or 4 . 0 , - 2 , 1*3. , 

e. You are then asked if you want to see a TECPLOT display. If you answer yes, follow 

the instructions in the next section. 

f You are asked to verify if the surface is correct. If not, you are prompted to enter data 
for the surface again. 

g. You are asked if you want to rotate or translate the surface. If yes, you are grven the 
following menu: 

SELECT ACTION 
T TRANSLATE SURFACE 

X ROTATE ABOUT X AXIS 

. Y ROTATE ABOUT Y AXIS 

Z ROTATE ABOUT Z AXIS 

h. Enter the menu choice, and then how you much you want the surface rotated or 

translated 5 . , 

i. Again, you are asked if you want to see a TECPLOT display. If yes, follow the 

instructions in the next section. 

j. At this point, you have the chance to undo the operation if it doesn't look right. 

k. Further translation and rotation is available at this point 

l. Finally, the main menu is displayed again. 

7. To modify an existing surface: 

a. Select 2. 

b. A list of surfaces is displayed. Choose the number of the surface you want to modify. 

c. The following menu is shown: 

SELECT AN ACTION BY NUMBER 

1 MODIFY DATA FOR THE SURFACE 

2 ROTATE OR TRANSLATE THE SURFACE 

3 DISPLAY GEOMETRY AND SURFACE DATA FOR THE SURFACE 

4 FINISHED MODIFICATION OF THIS SURFACE 

d. Choose the appropriate actions until the surface is to your satisfaction. At various 
times, you will be given the option to view the surface in TECPLOT. When finished 
with the surface, choose number 4. The main menu will be displayed. 

8. To delete a surface: 

Note: Be careful - this operation cannot be undone. 

2 L Select 3 . 

b. A list of surfaces is displayed. Choose the number of the surface you want to delete. 

c. The surface is deleted. The main menu will then be displayed. 

5Note that rotatio ns actually take place about an axis through a user-selected point parallel to the selected coordinate 
axis of rotation rather than the coordinate axis itself. This allows the user great flexibility when rotating surfaces. 
For example, consider a rectangle initially in the Y-Z plane and intersecting the X axis at xq as shown in figure 
2.2.3-l(a). If this rectangle is rotated by 90 degrees about the Z axis, with (0,0,0) as the rotation point, it is 
transformed to the position shown in figure 2.2.3- 1(b). However, if the rotation point is taken to be (xo,0,0), the 
rectangle is transformed to the position shown in figure 2.2.3-l(c). 
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9. To insert a new surface: 
si Select 4* 

b. A list of surfaces is displayed. Choose the number of the surface in front of which you 
would like the new surface to be inserted 

c. Enter information just as in step 6. 


10. To review surfaces: 

Select 5. 

b. A list of surfaces is displayed Choose the number of the surface you want to review. 
Zero reviews all surfaces. 

c. You are shown the properties of that surface, and then given the option to get a 

TECPLOT display of the structure. 

d. If you chose to review all surfaces, the program steps you through each surface, asking 
if you want to see a TECPLOT display at any point 

11. To copy a surface: . . , ... 

Note: This menu option is very useful for building structures containing surfaces which 

are identical except for positions. 

a. Select 6. „ 

b. A list of surfaces is displayed Choose the number of the surface you want to copy. 

c. A list of surfaces is displayed Choose the number of the surface after which you want 
the new surface to be inserted. Zero inserts the new surface before the first surface. 

d. You are then asked for a comment line for the new surface. 

e. At this point you have the option of modifying the surface, so that it can be rotated or 

translated. Follow the instructions in step 7. 

12. If you are editing a SHADOW file in order to create a SOLSHAD file, you will then be asked 
for the solar surface properties of all the surfaces that previously only had atomic oxygen 
properties. You will be given the opportunity to view the structure with TECPLOT, change 
the material name, and then enter the reflectivity and absorption coefficients. 

13. At this point, you are given the chance to review or modify any surfaces. This is particularly 
useful if you made a mistake in entering the solar properties. If you choose Y, then you are 
returned to the main menu. 


14. When finished choose 7 from the main menu. 

15. You then have the option of translating or rotating the entire structure. Enter the menu choice, 
and then how you much you want the surface rotated or translated. 

16. If you chose to rotate or translate, you then have the options of viewing the transformed 
object, and then undoing the translation or rotation if you desire. 

17. You may continue to translate or rotate the entire structure until it is where you want it to be. 
At this point you have the option once more to review all surfaces. 

18. The last section is for the input of various values needed in the calculation. The program 
displays the current values and asks if you want to modify them. If this is a new file, then the 
current values are generally zero. 

19. Enter the start year, day, and hour. First enter the year and day in the format YYYY DDD, 
where the day is the number of days from the beginning of the year. 

20. Enter the Greenwich Mean Time hour and fraction. 
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21. Enter the end year, day, and hour in a similar manner. 

22. Enter information on the following values (see sec. 2.2.1 for more information): 

Range of angular positions on orbit 
Ascending node longitude and range of precession 

23. Enter whether the event time and orbit precession are correlated or uncorrelated. 

24. Enter information on the following values (see sec. 2.2.1 for more information): 

Orbit inclination 
Earth albedo 
Earth radius 

Average circular orbit satellite altitude 

Maximum number of direct solar rays to scatter 

Maximum number of Earth-reflected solar rays to scatter 

Number of Sun positions after which Earth-reflected radiation will be calculated 

Number of randomly selected Sun positions 

25. Enter if you want the surfaces to be flagged if rays intersect with inactive sides during 
generation of the visibility matrix. 

26. Enter information on the following values (see sec. 2.2.1 for more information): 

The number of zenith and azimuth increments 
Random number generator seed 

27. Enter the name of the file to which the data will be written. 

28. Program is complete. Click on “Yes” to exit the window. Return to the TECPLOT window 
(Type <CntlxEscape>. Double-click on Tecplot.) Quit (Type: /fqq) 

Viewing a TECPLOT Display of the Surfaces 

1. Wait until the message “SWITCH TO TECPLOT NOW AND DISPLAY SURFACES” is 
written. Then type <CntlxEscape>. The box with available windows will appear. 
Double-click on Tecplot. 

Note: With Windows 3.1, you can use <AltxTab> to switch between MDDB and 
TECPLOT once the first switch has been established using <CntlxEscape>. 

2. You may get a message saying that there is insufficient memory for the display. If this occurs, 
click OK, and then double-click on the Tecplot icon at the bottom of the screen. (What this 
means is that when TECPLOT first appears after switching from another window, there will be 
some extraneous garbage on the screen. However, this will disappear as soon as the screen is 
refreshed.) 

3. Type <CntI>p p p<enter>. This will play a macro called p . mcr, which takes the file 
mddb .pit generated by MDDB and plots in on the screen. (If there is extraneous garbage 
on the screen, this step will eliminate it.) Each of the surfaces will be a different color, and 
arrows will appear that show which side of each surface is the active side. 

Note: If an error occurs while trying to run p . mcr, this may be due to the binary file format 
generated by the FORTRAN code not being compatible with the TECPLOT version in use. To 
fix this, quit TECPLOT, create a DOS shell, copy p2 .mcr to p .mcr, exit the DOS shell, and 
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start TECPLOT again. Now when you type <CntI>p p p<enter>, the macro will run the 
ASCH file through PREPLOT to generate the binary file. (p2 . mcr will take a little longer to 
run than p .mcr.) 

4. To rotate the object or identify surfaces, see the directions below. If the display appears 
incomplete, refer to the section on possible problems with TECPLOT. 

5 To return to MDDB and continue, type <CntlxEscape>, and double-click on 

MDDB - [Unit *]. 

Note: With Windows 3.1, it is possible to switch back and forth between windows by typing 
<AltxTab>. Refer to the Windows manual. 

Rotating the Structure in TECPLOT 

If the view that TECPLOT gives you is not exactly what you want, it is a simple matter to 
rotate it Type the following commands: 

/ Get to main menu, 

c Get to Contour menu 

v Get to view menu 

r Rotate 

x,y, or z Choose x, y, or z axis to rotate about 

Now use the mouse or the keyboard arrows to rotate the object <Esc> will bring you out 
of rotation mode. If the rotation has moved some parts of the object outside of the field of view, 
type v (for view menu) and then f (for fit). 

Type r to regenerate the picture. 

Note: Every that MDDB generates a new picture for TECPLOT, it reverts back to its original 
rotation view. To save a particular view, store it as a style sheet. Once you have the view 
you want, type /fsw and the filename (such as mddb . sty). To restore this view later on, 
type /fsr <filename>. 

Determining Surfaces in TECPLOT 

The TECPLOT display shows all the surfaces in different colors so that you can determine 
which surface on the screen corresponds to which surface in the SOLSHAD or SHADOW input 
file. There are two ways to make this determination: displaying a table of colors with the surfaces 
they correspond to, and individually examining surfaces. 

To generate a table, type the followi ng commands: 

1. <Cntl>p p table<enter>. (This plays a macro called table .mcr.) 

2. Wait until the message “Enter total number of surfaces” appears on the bottom of the screen. 
Press <enter>. 

3. Type the total number of surfaces. (This was displayed in MDDB just before transferring to 
TECPLOT.) Press <enter>. 

4. Wait until the message “Place the table” appears at the bottom of the screen. Use the mouse to 
put the table where you want if Click with the left button to place it. You can read the table by 
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matching up the color with the number just directly above it. (The one exception is surface #1, 
which is always at the very bottom of the table.) 


To individually examine surfaces, type t he following commands: 

/ Get to main menu. 

c Get to Contour menu 

p Get to probe menu 

e Examine 

Now use the mouse to place the crosshairs on the region of interest. Click with the left 
mouse button. A table will appear on the left side of the screen. The second value in the table 
corresponds to the surface number. (If you get a non-integer number, it is because you have 
clicked on an arrow instead of a real surface.) To get out of examine mode, hit <escape>. To get 
rid of the crosshairs, regenerate the screen (/or). 

Possible Problems 

A minimum number of grids are required for certain shapes: for example, you cannot create 
a disk with only 2x2 grids. The second grid value NN must be at least 3 for TECPLOT to 
approximate the disk by a triangle; use a value of 8 or more to get anything remotely circular. 

TECPLOT has difficulty with mesh units that come to a point, and it will generally not 
draw them. To avoid this problem, do not use values of zero for ratios in the trapezoid, disk, and 
cone geometries. Instead, use ratios that are very small (0.001) to approximate triangles (ratio of 
trapezoid short size to long side), disks with no holes in the center (ratio of inner radius to outer 
radius), and cones that come to a point (ratio of small radius to large radius). Other coordinates 
may haVe to be modified slightly if you are concerned about geometries overlapping. Similarly, 
spheres’ polar angular extent should range from 0.001 to 179.999 instead of 0 to 1 80 degrees. If 
you choose not to add this small factor of 0.001, SHADOW will still perform the calculations 
correctly, but the TECPLOT display may be incomplete, both during input file generation and 
during die final display of solar exposures. 

TECPLOT will occasionally have other quirks as well. Sometimes mesh lines are missing 
or are drawn where they shouldn't appear. In general, rotating the view to something slighdy 
different will correct these problems. 


2.3 INSTALLATION AND OPERATION OF THE MICROENVIRONMENT 
SOLAR EXPOSURE PROGRAM 

This section describes installation and operation of SOLS HAD on the CONVEX computer 
on which it was developed, and gives some suggestions for migrating SOLSHAD to other 
computers. 

2.3.1 Installation on a CONVEX Computer 

SOLSHAD was developed on a CONVEX C2 supercomputer running under CONVEX 
operating system release 9.0 and CONVEX version 4.2bsd unix. SOLSHAD is written in 
CONVEX FORTRAN. This section describes the installation and operation on a machine identical 
to that on which it was developed. It is suggested that a directory called sol micro and one 
subdirectory called solshad be created, and that all files be put into these directories. In what 
follows it is assumed that these directories have been created. 
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All files for SOLSHAD reside on directory solmicro and its subdirectory solmicro/solshad. 
Subdirectory solmicro/solshad contains FORTRAN source code for program SOLSHAD and its 
subroutines as well as a Makefile for building the SOLSHAD executable. These files are listed in 
tables 2.3.1-1 and 2.3.1-2, respectively. 

Note : Users of the microenvironment atomic oxygen exposure program SHADOW (ref. 1) will 
note that SOLSHAD and SHADOW share a number of subroutine names in common. Although 
their names are the same, many of these subroutines are different and cannot be interchanged 
between the two codes. 


Table 2.3.1-1 

Files in the Directory solmicro 

run.solshad* solshad.in_sample solshad.tp7_sample 

solshad/ solshad.out_sample solshad.x 


Table 2.3.1-2 

Files in the Directory solmicro/solshad 


Makefile 

cwb2int.f 

hiterth.f 

rayg3.f 

switch.f 

autovis.f 

cwbint.f 

les.f 

rayg4.f 

unitized 

ceshint.f 

cwbsc.f 

multi.f 

reflectd 

xrotd 

cload.f 

dateda.f 

packer.f 

rglerr.f 

yrotd 

connec.f 

diffused 

pmote.f 

solshad.f 

zrot.f 

cwb.f 

dot.f 

raygl.f 

stored 


cwb2.f 

doy.f 

rayg2.f 

subsold 



Follow the procedure below for loading SOLSHAD and preparing it to execute. In this 

procedure, bold text indicates commands the installer should type. Press return after each 

command. Remember that the unix operating system is case sensitive. 

1 . C reate and load the directory solmicro with the files given in table 2.3. 1 - 1 and subdirectory 
solmicro/solshad with the files in table 2.3. 1-2. 

2. Get to subdirectory solmicro/solshad by typing cd /pathname/solmicro/solshad, where 
pathname is the appropriate path. 

3. Type make. This causes the file Makefile to be run, compiling all FORTRAN routines in 
the subdirectory and creating an executable called solshad.x. 

4. Type cp solshad.x .. to copy the executable to solmicro. 

5. Get to directory solmicro by typing cd .. 

6 . Edit file run.solshad to change the lines 

solshad.x < solshad.in_sample >! solshad.out_sample 
to 

solshad.x < solshad.in_sample >! test.out 
and 

mv TAPE7 solshad.tp7_sample 
to 

mv TAPE7 solshad.tp7_test 
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7. Change the path in the second line to the proper path. Save and exit. 

8. Ensure that run.solshad has execute permission by typing chmod +x run.solshad. 

9. Run the sample test case by typing run.solshad. The sample test case takes about 3.5 
minutes CPU time. When the run is complete, compare file solshad.out_sample with 
test.out and solshad.tp7_sample with solshad.tp7_test. These files should give the same 
results except for minor differences in CPU times. 


2.3.2 Running the Microenvironment Solar Exposure Program 

Before running SOLS HAD, build its input file describing your geometry and operating 
conditions using MDDB (per sec. 2.2). Edit the run.solshad file in the solmicro directory to direct 
your program input file to standard input and to direct standard output to your output file. You 
may also wish to direct TAPE7, the file containing information for displaying the solar exposure in 
TECPLOT, to a different file name. Every time SOLS HAD is run, any existing version of TAPE7 
is overwritten. SOLS HAD may be run interactively by typing run.solshad or in batch mode by 
typing qsub -q v run.solshad where v is the queue to which the job is to be sent. Because 
SOLSHAD may take significant execution time, running in batch mode is recommended. 

As written, SOLSHAD is set to allow up to 100 surfaces and 5,000 nodes. These 
dimensions are set by PARAMETER statements in the FORTRAN code. Table 2.3. 1-3 
summarizes these parameters and what they affect Users who do not run cases which need these 
large dimensions may find that SOLSHAD runs more quickly if they are reduced. 

Table 2.3.1-3 
Parameters in SOLSHAD 



Value 


EORTRAN " 

i;i Mi'll Si H M 



roAmm:! 

Routines Affected 

usage 

Block 



MAXPS 

100 

//, /SURFC7, 

/SURF/, 

/SEE/ 

AUTOVIS, CESHINT, 
CLOAD, CONNEC, 
MULTI, RAYG1, 
RAYG2, RAYG3, 
RAYG4, SOLSHAD 

Maximum number 
of surfaces. 

MAXNOD 

5000 

//,/FLX/, 

/RNODE/, 

/RTR/ 

AUTOVIS, CESHINT, 
CLOAD, CONNEC, 
CWB2INT, MULTI, 
PACKER, RAYG1, 
RAYG2, RAYG3, 
RAYG4, SOLSHAD 

Maximum number 
of nodes. 

MAX© 

3*MAXNOD 

1 

// 

AUTOVIS, CESHINT, 
CONNEC, MULTI, 
RAYG1, RAYG2, 
RAYG3, RAYG4, 
SOLSHAD 

Dimension of © 
array. 


• All FORTRAN routines are in files of the same name in lower case with a .f extension. 

• All files are are in directory solmicro/sol shad. 
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2.3.3 Some Information for Installation on Other Computers 

This section provides information about SOLSHAD which will be useful if SOLSHAD is 
to be converted to run on computers other than the CONVEX described in section 2.3.1. 

SOLSHAD is written in mostly ANSI standard FORTRAN and generally does not call 
functions specific to any particular operating system or FORTRAN compiler. However, the 
following exceptions arc noted. 

Calls are made to CONVEX lveclib8 CPU elapsed time function CPUTIME in program 
SOLSHAD. This function must be changed to the proper function name for other systems, or, 
references to it may be deleted without significant effect on program execution. CPUTIME is used 
to output elapsed CPU time at various stages of SOLSHAD. 

RAN is a CONVEX lveclib8 VAX-like random number generator. The appropriate 
random number generator call must be substituted in subroutines DIFFUSE, MULTI, RAYG4, 
and in program SOLSHAD. Also, the random number seed initialization must be made proper for 
the random number generator used. RAN uses a large integer number as a seed. RAN generates 
uniformly distributed random numbers between 0 and 1. 

Numerous seven-character variable names are used in SOLSHAD. It is believed that most 
can be truncated to the ANSI standard six characters without causing duplicate variable names. 
However, care should be taken before doing this to ensure that there are indeed no duplicate 
variable names generated. 

SOLSHAD is compiled to run using CONVEX double-precision word length (8 bytes or 
64 bits) for all real variables. This word length is needed to maintain numerical accuracy. Real 
variables in FORTRAN routines in solmicro/solshad are not declared double precision explicitly, 
but are converted to double precision with a FORTRAN compiler option. Users converting 
SOLSHAD to run on other computers should be aware that single-precision word length (4 bytes 
or 32 bits) real variables may not give the same numerical accuracy as double-precision word 
lengths. 

The following CONVEX FORTRAN compiler options are either recommended or 
mandatory when compiling SOLSHAD: 

-or all Provides full optimization report (recommended). 

-02 Optimize code with vectorization (recommended). 

-pd8 All default integer, logical, real, and double-precision values occupy 8 bytes 

of storage (mandatory). 

-rl Perform loop replication optimizations (loop unrolling, dynamic code 

selection) when profitable to do so (recommended). 

-72 Process only the first 72 characters of each line. Lines containing tab 

characters will not compile properly with this option (mandatory). 

2.4 OUTPUT: 3D PLOT GENERATION 

2.4.1 Plot Generation 

SOLSHAD generates an ASCII file that can be read by TECPLOT to display the structure 
with colors on its surfaces that correspond to the different levels of CESH solar exposure. This 
file is named TAPE7, unless otherwise renamed in the batch file. First, the ASCII file must be 
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converted into a binary file for input into TECPLOT, using a program called PREPLOT. Then a 
macro can be run to display the structure. Finally, TECPLOT commands can be used to rotate, 
add text, print, extract data, etc. 

For complete instructions on what can be done in TECPLOT, read the software manual. 
This section contains instructions on how to do the most common tasks for viewing the 
SOLSHAD output using a PC. To get started, follow these instructions: 

1. Change directory to where TECPLOT is located. (Generally, cd \tecplot.) 

2. Type preplot <fllename.xxx>, where <filename.xxx> is the TECPLOT file generated by 
SOLSHAD. PREPLOT converts the file into binary, and names die new file <filename.plt>. 

3. Type tecplot to start up TECPLOT. 

4. Load the data file. (Commands: f d r. Type in the file name with .pit extension.) 

5. Run the macro shadow. mcr. To do this type <cntl>p p shadow. The structure will 
appear, although it may be difficult to see various parts of it due to the positioning of the x, y 
and z axis. 

2.4.2 Plot Manipulation 

Rotating the Structure 

Type the following commands: 

/ Get to main menu, 

c Get to Contour menu, 

v Get to View menu 

r Rotate. 

x,y, or z Choose x, y, or z axis to rotate about. 

Now use the mouse or the keyboard arrows to rotate the object <Esc> will return bring 
you out of rotation mode. If the rotation has moved some parts of the object outside of the field of 
view, type v (for View menu) and then f (for Fit). 

Type r to regenerate the picture. 

Adjusting the Colors and Adding a Color Table 

TECPLOT will automatically choose a color scale, based on the highest and lowest values 
of the CESH. In general, this scaling is good, but sometimes a small number of nodes will have 
extreme values that will throw off the scale. To adjust the color scale, type the following 
commands: 

/ Get to main menu, 

c Get to Contour menu, 

c Get to Contour Value menu. 

I Levels, 

n New Range, 

n New Minimum and Maximum. 
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Choose a number of levels. (15 is a good value.) Enter the minimum. Enter the 
maximum. Type <cntl>r to regenerate the picture. Redo the steps from “New Range” to adjust 
until you get a good set of contours. 

To add a table that shows which colors correspond to which levels, type the following: 

/ Get to main menu. 

c Get to Contour menu. 

c Get to Contour Value menu. 

t Table. 

In this menu, you may adjust the text fonts and colors, the spacing in the table, and 
whether you have a bar outline. Finally, type p to place, and use the mouse to put the table where 
you want it to be. Click the mouse to finalize die table location. 

Plotting Specific Ray Groups 

The total CESH on a surface can be divided into primary CESH, specularly-reflected 
CESH, diffusely-reflected CESH, and absorbed CESH. (These groups can then be divided into 
direct solar CESH and Earth-reflected solar CESH.) It is possible to view the structure looking at 
only exposure from each of these groups. Table 2.4.2- 1 shows the variable numbers and names 
that correspond to the various ray groups. 

Table 2.4.2-1 

TECPLOT Variable Numbers and Names 


Variable number Variable name Value 


1 

NODE 

2 

PRIM 

3 

SPEC 

4 

DIFF 

5 

ABSOR 

6 

TOTAL 

7 

PRIMA 

8 

SPECA 

9 

DIFFA 

10 

ABSA 

11 

TOTAL 

12 

GRTOT 

13 

X 

14 

Y 

15 

Z 


Node index 

Primary direct solar CESH 

Specularly-reflected direct solar CESH 

Diffusely-reflected direct solar CESH 

Absorbed directsolar CESH 

Total direct solar CESH 

Primary Earth-reflected solar CESH 

Specularly-reflected Earth-reflected solar CESH 

Diffu sely-reflected Earth-reflected solar CESH 

Absorbed Earth-reflected solar CESH 

Total Earth-reflected solar CESH 

Total of direct and Earth-reflected solar CESH 

x-coordinate 

y-coordinate 

z-coordinate 


Type the following commands: 

/ Get to main menu. 

c Get to Contour menu. 

c Get to Contour Value menu. 

v Get to tiie Variable menu 

Use arrows to choose your CESH group. Press enter 

<esc> Get to Contour menu 

r Regenerate the image 
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Style Sheets 


A style sheet contains all the information as to how a data file is displayed It can be used 
on the data that it was generated with or on new data that is similar in format. To save or restore a 
style sheet, type the following: 

/ Get to main menu, 

f Get to File menu, 

s Get to StyleSheet menu, 

w or r Write or Read 

filename Filename that ends with .sty. 

Printing 

Type the following commands: 

/ Get to main menu, 

f Get to File menu, 

p Get to Print menu. 

In this menu, you may adjust the paper size and orientation, the device configurations, the 
printer format, number of copies, the file path for hardcopy files, and where to route the output. 
When done, type g for GoPrint Now you can position and scale the picture. Again, type g for 
GoPrint when done, and the print will be made. 

Probing Data 

TECPLOT allows you to pick a spot on the surface and get the exact exposure value there. 
Type the following commands: 

/ Get to main menu, 

c Get to Contour menu, 

p Get to Probe menu, 

e Examine. 

Use the mouse to place the cursor where you want the exposure value. Click with the left 
button. A table appears on the left side of the screen. You can read off the values next to the 
variable numbers. (See table 2.4.2- 1 for an explanation of these numbers.) You can continue 
clicking on new points. Press esc or right mouse button when done. 

Extracting Data 

TECPLOT allows you to draw a path across your surface, extract the CESH values, and 
generate a 2D plot of CESH versus distance. Type the following co mmands - 

/ Get to main menu, 

c Get to Contour menu, 

p Get to Probe menu, 

x Extract data. 

Use the mouse to place the cursor at the beginning of the path. Click with the left button. 
Move the mouse to the end of the first line segment. Click with the left button again. Continue 
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until you get to the end of the path. Then click the right mouse button. Enter y to accept the 
polyline, or n to redo it. Enter the number of points for your x-y curve. Enter a filename and a 
header. 

A new file has been created with this data. You can look at it by loading the new file. 
(Commands: / f d r. This new file will be automatically inserted as the default name.) Plot it as 
an X-Y plot. (Commands: / x). You will need to set the variables (x should set to the X, Y, or Z 
variable, and y should be set to an exposure value, such as GRTOT. Command: / x o v) and refit 
the graph to the axis (Command: / x v f)- 
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3.0 TECHNICAL FEATURES OF THE COMPUTER MODEL 


3.1 PROGRAM FLOWS 

3.1.1 Program Flow for the Microenvironment Solar Exposure Model 

Figure 3.1. 1-1 shows the SOLSHAD block diagram. The names of major subroutines in 
the block diagram are given in parenthesis. Figure 3. 1.1 -2 is a subroutine tree for SOLSHAD and 
shows all subroutines and functions and their calling hierarchy. This section describes the general 
order of program execution; details of the algorithms used are given in section 3.2.1. 

SOLSHAD starts execution by echoing its input file to output The input file is then 
rewound. Subroutine RAYG1 reads the geometry definition and surface properties for all surfaces 
and sets up the grid of nodes for each surface and other quantities needed for ray tracing. 
Calculation control parameters, including event start and end dates, satellite orbit description, Earth 
albedo, the number of Sun and orbit positions for which to trace rays, and whether multiple 
scattering is required, are read. 

All surfaces have previously been divided into a grid of nodes by subroutine RAYG1. 
SOLSHAD defines a set of points at which solar exposure in CESH will be calculated so that 
nearly the entire surface is covered. This grid for solar exposure has points at one comer (call it 
the lower left) of each node plus points in the lower right comer of the right-most column of 
nodes, at the upper left of the topmost row of nodes, and in the upper right comer of the upper 
right node. Each point is set in 1 % of the node size from the nearest edge of the node. This gives 
nearly full coverage of the surface and avoids any ambiguity as to which node the point belongs to. 
The coordinates of each grid point in Cartesian (x, y, z) space are calculated by subroutine 
RAYG2. Subroutine RAYG3 calculates the surface normal vector and two orthogonal tangent 
vectors. These data are saved for repeated use in subroutine AUTO VIS. 

Subroutine AUTOVIS calculates a visibility matrix for the surfaces defining the geometry 
of the object The visibility matrix indicates which surfaces can be shielded by themselves or other 
surfaces, or can reflect solar radiation to other surfaces or to themselves. The visibility matrix is 
analyzed to determine which, if any, surfaces have solar radiation reflected to them by any surface. 
This information is used to speed up program execution in two ways. First, when surfaces are 
being checked to see if a ray to or from them will be blocked by any surface, only those surfaces 
which can see the first surface need be checked. Other surfaces can be ignored. Second, if a 
surface cannot have solar exposure reflected on to it, there is no need to ray trace reflected solar 
exposure from that surface. 

The visibility matrix is calculated by ray tracing from each point on every surface to every 
other point and ray tracing along the surface normal from the active side of the surface at each 
point The process is sped up by noting that plane surfaces (trapezoids and disks) and the positive 
sides of cones, cylinders, and spheres can never reflect onto themselves. During this ray tracing, 
rays which hit the inactive sides of surfaces are flagged. Such rays often indicate that the object 
geometry is not properly constructed and that surfaces may be exposed to solar radiation which 
were not intended to be exposed. 

After the visibility matrix has been calculated, data for the centers of all the nodes in the 
structure are calculated and saved for repeated use during CESH calculations later in the program. 
These data include the coordinates of the center points and the surface normals and tangent vectors 
at the points. Also at this time the average area of a node is calculated. The average area is used to 
determine how many rays of CESH to scatter from each node. 
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SOLSHAD BLOCK DIAGRAM 



START 


Read geometry input and 
precalculate ray trace parameters 
(RAYG1 ) 


I 


Read calculation control parameters 


Precalculate grid point coordinates, 
surface normal, and tangent vectors 
(RAYG2, RAYG3) 


Form visibility matrix 
(AUTOVIS) 


Save data for node centers 
and get average node area 
(RAYG2, RAYG3) 


Loop over Sun positions § 

Choose random time, orbit position, and Sun position 
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Generate vector from structure to Sun 
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Figure 3. 1.1-1. SOLSHAD Block Diagram (Sheet 1 of 4). 
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(CESHINT) 
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Figure 3 . 1 . 1 - 1 . SOLSHAD Block Diagram (Sheet 4 of 4). 
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SOLSHAD Subroutine Tree 


SOLSHAD- 


PRNOTE 

RAY61 - 

IDOY 

DATEDA 

RAYG2 

RAYG3 

AUTOVIS 


-DOY 

- SUBSOL 
-XROT 

- ZROT 

- ACOSD 

- ATAN2D 

- CLOAD 

- DOT 

- RAYG4 - 
-MULTI - 


— HITERTH 
— CESHINT 


l— CONNEC 


RG 1 ERR 


CLOAD 

STORE 

SWITCH 

UNITIZE 

DOT 

RAYG4 - 


RAYG3 


RAYG3 

1 — SWITCH 

— DOT 
— RAYG2 
— RAYG3 
— REFLECT 
— DIFFUSE 
• — RAYG4 — RAYG3 

T PACKER 

CWB2 1— CWBC 

I— CWB 

CWB2INT r CWBINT 

LES L CWB 


Figure 3. 1.1-2. SOLSHAD Subroutine Tree. 


The portion of the surface of the Earth seen by the satellite is divided into a grid of points 
spaced equally in zenith and azimuth angles relative to a set of spherical coordinates centered at the 
satellite and with polar axis directed from the satellite to the center of the Earth. The sines and 
cosines of the angular coordinates of these points are stored for later used during ray tracing of 
Earth-reflected solar radiation to the satellite surfaces. 

At this point, SOLSHAD is ready to start calculating the solar exposure to the points of the 
surfaces of the object. An allowed random position along the satellite orbit is chosen. A random 
day of year and time of day are fed to subroutine SUBSOL, which returns the subsolar latitude and 
longitude on Earth; that is, the latitude and longitude at which the Sun is directly overhead. The 
Sun position coordinates are translated into the coordinates used by the object geometry and the 
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vector from the satellite to the Sun is saved. A check is made to see if the Earth is between the Sun 
and the satellite. If so, no solar exposure to the satellite occurs and new Sun and satellite positions 
are chosen. By choosing a large number of Sun and satellite positions, a statistically valid 
approximation of the solar exposure to the satellite is generated. 

A Sun and satellite position having been chosen that allows solar exposure to the satellite 
surfaces, each node on the object is considered in order. The procedure for every node on a 
surface after the first node on each surface is identical. For all direct solar illumination 
calculations, the center of the node is used as the proxy for the entire node. 

For the first node on a surface, flags are set indicating that ray tracing from the point at the 
center of the node to the Sun must be done to determine if the point is exposed to direct solar 
illumination. The previously calculated coordinates, normal vector, and tangent vectors of the 
surface at the point are retrieved. A check is made if the Sun is shining ot the active side of the 
surface. If not and the surface is a trapezoid or disk, the Sun cannot shine on any point on the 
surface, and flags are set to indicate that the surface is not exposed to solar illumination. In this 
case, no ray tracing is necessary for the surface. The rest of the procedure for determining solar 
exposure is identical for all nodes on the object 

If the node could perhaps be exposed to solar radiation, subroutine RAYG4 is called to 
determine if any surface blocks exposure. If the ray to the Sun is unblocked, direct solar exposure 
is accumulated. If selected, subroutine MULTI performs multiple scattering of the ray. The 
number of rays to Monte Carlo scatter is determined as the dot product of the surface normal at the 
point and the ray to the Sun times the number of rays per average node area times the ratio of the 
node area to the average node area. This ensures that each surface is illuminated with the same 
areal density of CESH rays. The scattering procedure is the same for each scattered ray. 

For each scattered ray, subroutine MULTI reverses the direction of the ray so that it is 
directed toward the center of the node rather than to the Sun. A check is made to see if the ray is 
parallel to the surface or very nearly so. If it is, the ray is given no further consideration; the ray 
tracing routine has difficulty handling such rays properly and, in any event, they contribute only a 
very small solar exposure. Rays not parallel to the surface are Monte Carlo scattered by the 
following procedure. 

A ray of solar CESH striking a surface may undergo one of three fates: specular reflection, 
diffuse reflection, or absorption. In the event of absorption, propagation of the ray stops and the 
exposure due to the absorption is accumulated. The fate of a ray is determined by a Monte Carlo 
(random) selection in which the probability of selection is proportional to the surface property 
value of the fate. The surface properties sum to one, so one of the three fates is guaranteed to be 
selected. 

If reflection is selected, a random point on the node is chosen and its coordinates, surface 
normal, and tangent vectors are calculated by RAYG2 and RAYG3. For nonplanar surfaces, the 
surface normal is a function of position on the surface and the tangent vectors may also be 
functions of position. Choosing random positions on a node for reflecting CESH ensures that 
surfaces which are exposed to reflected CESH receive the proper exposure independent of their 
area or orientation with respect to the ray direction. 

If specular reflection is selected, the normal component of the ray is reversed to give the 
reflected ray. Diffuse reflection has a probability of scattering uniform in cosine of zenith angle 
measured from the surface normal and uniform in azimuth angle. These two angles are chosen 
randomly and the ray is directed along them. Subroutine RAYG4 propagates the reflected ray and 
determines whether it strikes the active side of a surface or not. If the ray does not strike the active 
side of a surface, consideration of that ray is terminated. If the ray does intersect the active side of 
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a surface, subroutine RAYG4 determines the location (node) on the surface and subroutine 
RAYG3 determines the surface normal and tangent vectors at the intersection. The specular or 
diffuse solar exposure at this point is assigned to the node and accumulated there. The ray itself is 
Monte Carlo scattered from the point of intersection on the node and propagated or terminated in 
the same manner as for the initial scattering. This process continues until the ray is removed from 
consideration or until a maximum of 100 scatters has occurred. 

After primary and direct reflected solar exposure have been calculated for all nodes. Earth- 
reflected solar radiation to the object is calculated if the Earth albedo is nonzero and calculation of 
Earth-reflected radiation is selected. The portion of the Earth viewed by the satellite is divided into 
a grid of points. The calculation of Earth-reflected solar exposure to the satellite from each point 
takes approximately as long as the total direct solar exposure calculation. Because of this, the user 
is given the option of calculating Earth-reflected solar exposure after every NALB^h direct solar 
calculation, where NALB may range from 1 to the total number of Sun positions. In the former 
case, one Earth-reflected solar exposure calculation is done per Sun position; in the latter, only one 
Earth-reflected solar exposure calculation is performed. If NALB is greater than the total number 
of Sun positions, no Earth-reflected solar exposure is calculated. The procedure for calculating 
Earth-reflected solar exposure from each grid point on Earth is the same. 

A unit vector ray from the satellite to the grid point on Earth is formed and saved. The ray 
is propagated to Earth and the outward normal vector from the intersection point with the surface of 
the Earth is calculated. The dot product between this normal and the ray to the Sun determines 
whether the point on Earth is illuminated by the Sun. If it is not illuminated, no Earth-reflected 
solar exposure calculation need be done. Otherwise the solar radiation from the Earth grid point is 
deterministically diffusely reflected toward the satellite. This diffusely reflected solar radiation is 
the incident solar exposure multiplied by the albedo times the projected area of the Earth grid 
element times the solid angle of the Earth grid element as seen by the satellite. Exposure to satellite 
surfaces using this Earth-reflected solar radiation is calculated in exactly the same manner as direct 
solar exposure to satellite surfaces. Direct Earth-reflected solar exposure, specular and diffuse 
Earth-reflected (from one satellite surface onto another) solar exposure, and absorption of Earth- 
reflected radiation are accumulated. 

When all solar exposure calculations have been completed, subroutine CESHINT 
interpolates the CESH (which to this point have assumed to be calculated at the center of the nodes) 
to the comers of the nodes using Stineman's consistently well-behaved method of interpolation 
(ref. 4). This gives full coverage of the surfaces to ensure proper looking plots of the CESH in 
TECPLOT. 

SOLSHAD generates a table of direct primary, specular-reflected, diffuse-reflected, 
absorbed, and the total of primary and specular- and diffuse-reflected solar exposure at each point 
on the structure, as well as the same for Earth-reflected solar exposure. This table is written to 
standard output Execution of SOLSHAD concludes by calling subroutine CONNEC to generate a 
text file suitable for plotting by TECPLOT after processing by the TECPLOT PREPLOT module. 
This file contains the various solar exposures in CESH at each point and a connectivity matrix 
which tells TECPLOT how to connect the points making up the structure. 

3.1.2 Program Flow for MDDB Version 2.0 

Figure 3. 1.2-1 shows the MDDB block diagram. The names of major subroutines in the 
block diagram are given in parenthesis. Figure 3. 1.2-2 is a subroutine tree for MDDB and shows 
all subroutines and functions and their calling hierarchy. The majority of MDDB execution 
involves user responses to prompts printed to screen by MDDB. Major computational routines to 
specify surfaces point by point from the geometry definition such as CONNECM, RAYG1M, 
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RAYG2M, and RAYG3M are stripped down versions of the subroutines of similar name in 
SOLS HAD. Consequently, detailed description of the algorithms used in MDDB is unnecessary. 


MDDB Block Diagram 
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Figure 3. 1.2-1 . MDDB Block Diagram (Sheet 2 of 3). 
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Figure 3. 1 .2-1. MDDB Block Diagram (Sheet 3 of 3). 
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MDDB Version 2.0 Subroutine Tree 
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Figure 3. 1.2-2. MDDB Subroutine Tree (Sheet 1 of 2). 
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Figure 3. 1.2-2. MDDB Subroutine Tree (Sheet 2 of 2). 


MDDB starts execution by asking the user whether he or she wishes to create an input file 
for SHADOW Version 1 .X or Version 2 X or SOLS HAD and whether to create a new input file or 
to edit an existing one. If the user chooses to edit an existing input file, MDDB asks the user to 
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specify the program to which the file to be edited applies and asks for the filename. The data from 
the file are read in. At this point the user has access to edit all of the surfaces in the object In 
addition, the user may perform any of the file building and editing operations available when a new 
input file is being created. These operations are described below. 

If a new input file is chosen, the user is prompted to enter the geometry and surface 
property definition for the first surface. When the data have been entered, the user is given the 
opportunity to display the surface using TECPLOT. When display mode is selected, the geometry 
specifications are written to a scratch file on unit 11. This file is read by subroutine RAYG1M 
which does preliminary setup and interpretation on the geometry inputs. Subroutine CONNECM 
calls subroutine RAYG2M to locate the Cartesian coordinates of each point on the surface of the 
structure. These points are written to a binary file directly readable by TECPLOT and to an ASCII 
file which may be used by PREPLOT. Outward normal arrows from the active side of surfaces are 
also generated and written to these files. Last, the finite element connectivity matrix specifying 
how points on surfaces are to be connected are generated and written to the files. The TECPLOT 
display shows all surfaces comprising the structure with direction arrows pointed outward from the 
active side of each surface. 

The user may choose to rotate or translate the surface and look at it again with TECPLOT. 
If the user does not like the operation just performed, he or she may undo it Surface rotations take 
place about rotation axes parallel to the Cartesian x, y, and z axes. The user specifies the position 
of these axes. This gives greater freedom to move a surface than a rotation about one of the 
coordinate axes. 

When the first surface specification is complete or if the user is editing an existing input 
file, a menu of actions is presented. Actions include adding a new surface after the last entered 
surface or inserting a new one between existing surfaces, deleting a surface, reviewing surface 
(listing them with an optional TECPLOT display of the structure), modifying an existing surface, 
or copying an existing surface. The user repeats these actions until he or she is satisfied with the 
geometry and surface properties of the structure. 

After the user has finished building the geometry for the structure, MDDB checks to see if 
surface properties are correct for the type of input file being built For example, if the user is 
building a SOLSHAD input file by editing an existing SHADOW input file, the surface properties 
must be changed even though the object geometry has not been changed. 

After the user has finished building the geometry for the structure, the entire structure may 
be rotated or translated in the same manner as described for individual surfaces. 

For SHADOW Version 1.X input files only, MDDB prompts the user to enter or edit the 
visibility file, which specifies which surfaces can see themselves or other surfaces. The user has 
similar options to those available during geometry building: he or she may display the structure, 
review the visibility matrix, modify it, or build it from scratch. 

Last, the user enters new parameters or modifies existing parameters which control the 
atomic oxygen flux calculations for SHADOW or the CESH calculations for SOLSHAD. MDDB 
finishes execution by writing the completed input file to disk. 

3.1.3 TECPLOT Macros 

There are three macros to generate plots in TECPLOT. This section briefly describes what 
these macros do. 
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File: p.mcr 

1. Load the file MDDB.PLT. 

2. Turn off graphics. 

3. Create a mesh plot in 3D. Assign the Z axis to the last variable. 

4. Assign the x, y, and z axes to the last three variables. 

5. Choose to make a hidden line plot 

6. Export the mesh so that it will show up in the contour plot. 

7. Create a contour plot using the surface number (second variable) as the contour variable. 

8. Set zones as filled. 

9. Rotate it 30 deg in the x, y, and z directions to give it better perspective. 

10. Make sure the rotated plot fits in the screen. 

1 1. Activate the graphics and generate the plot 

Hatg; If p2.mcr is used, then the first step is to run PREPLOT with the file MDDB.PRE to 
generate MDDB.PLT. Otherwise, it’s the same as p.mcr. 

File: table.mcr 

1. Choose new contour levels: maximum, minimum, and delta. 

2. Set minimum level to 1. 

3. Put the message "Enter total number of surfaces" at the bottom of the screen. 

4. Get user input for maximum level. 

5. Set delta level to 1. 

6. Place the table. 

7. Get user input on where to place table, while displaying the message "Place the table." 

8. Regenerate plot. 

File: shadow.mcr 

1. Turn off graphics. 

2. Create a mesh plot in 3D. Assign the Z axis to the last variable. 

3. Assign the X, Y, and Z axis to die last three variables. 

4. Choose to make a hidden line plot. 

5. Export the mesh so that it will show up in the contour plot. 

6. Create a contour plot using the total solar exposure CESH (fourth variable from the end) as the 
contour variable. 

8. Set zones as filled. 

9. Activate the graphics and generate the plot 

3.2 ALGORITHMS IN THE MICROENVIRONMENT SOLAR EXPOSURE 
PROGRAM 

3.2.1 SOLSHAD Algorithms 

This section describes the algorithms used in SOLSHAD. Each algorithm is described 
independently of the others. Refer to section 3.1.1 for information on where the algorithms are 
used in SOLSHAD. 

Sun Position Determination. A Sun position is chosen by randomly choosing a day 
of year and time of day as GMT hour and fraction. These data are input to subroutine SUBSOL, 
which was adapted from the LOWTRAN7 code (ref. 5). SUBSOL produces an approximate 
(within 1 degree) subsolar latitude and longitude. The equation of time (minutes) and declination 
angle (degrees) for the day of year (small year-to-year variations of these quantities are ignored) are 
solved by interpolation from 25-point lookup tables. The declination angle is the same as the 
subsolar latitude SOLAT. The subsolar east longitude (degrees) is given by 
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( 3 . 2 . 1 . 1 ) 


SOLONG = 180 - 15 (t + Et) 

where t is the GMT hour and fraction and E t is the equation of time solution for the day of year. If 
needed, 360 degrees is added to the longitude to keep it in the range 0 to 360 degrees east. 

Orbit Position Determination. All calculations in SOLSHAD assume a circular 
satellite orbit around a spherical Earth. Such an orbit can be specified by its altitude HSAT, 

i nclinati on a, and east longitude of ascending node (fig. 3.2. 1-1); that is, the longitude at which 
die orbit crosses the equator in a northerly direction. The position of the satellite on the orbit can 
be specified by the angle P between a line from the center of the Earth to the ascending nodal 
crossing and a line from the center of the Earth to the satellite position. In SOLSHAD the satellite 
position on orbit is restricted to be between two limiting angles, ORBPOS(l) and ORBPOS(2). 

The range of these angles is arbitrary. A whole orbit is specified by setting them to 0 and 360 
degrees. Part of an orbit may be specified by restricting their range. The instantaneous position 

angle (J is selected by uniform random choice of angle between ORBPOS(l) and ORBPOS(2). 

The satellite orbit plane processes with the precession rate dependent on orbit geometry. 
This precession ma y be accounted for by specifying the longitude of ascending node at the event 
start PRECES(l) and the range of precession during the event PRECES(2). The processed 
position of the orbit plane during an event may be generated in either of two ways. The processed 
position of the orbit plane is chosen randomly from the range PRECES(l) to PRECES(l) + 
PRECES(2) using a uniform distribution. This choice ensures a statistically uniform coverage of 
all possible processed positions of the orbit This is the preferred method for lengthy missions 
during which the processed orbit position is not well known. 

Alternatively, the orbit plane may be chosen to process at a constant rate from its initial 
position at event start to its final position at the event end. To do this, the orbit precession is forced 
to be proportional to mission time. This choice is appropriate for brief missions during which the 
precession of the orbit plane with time is well known. 

At this point, the orbit plane precession is in terms of Earth east longitude at event start; 
however, in the time t days since event start the Earth has rotated to the east 360 degrees/day times 
t days. This Earth rotation during time t is subtracted from the longitude at event start to give 
RANPRE, the Earth longitude at time t. 

Transformation of Sun Position to Satellite Coordinates. The satellite 
coordinate system is attached to the satellite and is defined for satellites which maintain a constant 
orientation with respect to an axis radially outward from the center of the Earth through the satellite 
(fig. 3.2. 1-2). In this coordinate system the origin is at the satellite position. The X axis is 
defined to be radially outward from the center of the Earth through the satellite. The Z axis is 
defined to be perpendicular to the X axis and in the satellite heading direction on orbit. The Y axis 
is defined so that (X,Y,Z) form an orthogonal right hand coordinate system. In this coordinate 
system the center of the Earth is at (-RE-HSAT, 0, 0). 


55 


N 



a = Orbit inclination 

p = Orbit position angle 

HSAT = Satellite altitude above Earth 
N = Longitude of ascending node 
O = Center of Earth 

P = Satellite position 

RE = Radius of Earth 


Figure 3.2.1-1. Satellite Orbit Geometry. 
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Figure 3.2.1-2. Coordinate System Definition. 

The Sun position vector, which is defined as a unit vector from the center of the Earth 
through the subsolar latitude and longitude (that is, to the Sun), must be translated to be defined in 
terms of the satellite coordinate system. 

First, consider a coordinate system with its origin at the center of the Earth, the Z' axis 
through the North Pole, the X axis through the satellite orbit ascending node longitude, and the Y 1 
axis such that (X\ Y 1 , Z') form an orthogonal right hand coordinate system. In this coordinate set 
the subsolar latitude is unchanged, but the subsolar longitude is modified by subtracting off the 
ascending node longitude RANPRE, so 

SOLAT = SOLAT (3 .2. 1 .2) 

SOLONG’ = SOLONG - RANPRE (3.2. 1 .3) 

where, if needed, 360 degrees is added to SOLONG' to put it in the range 0 to 360 degrees. Then, 
the unit vector to the Sun in the ' coordinate system is (bold print indicates unit vectors; all vectors 
described in this section are unit vectors) 

S' = cos(SOLAT) cos(SOLONG') i’ +cos(SOLAT) sin(SOLONG') j' +sin(SOLAT) k’ 

(3.2. 1.4) 
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or 


S' = S’x i’ + S'y j' + S'z k* (3.2. 1.5) 

where i\ j', and k' are unit vectors in the X\ Y\ and Z' directions, respectively. 

Next, S' must be rotated about the X' axis by the negative of the orbit inclination, 
-ORBING, so that the rotated Y' axis is in the orbit plane and the rotated Z' axis is perpendicular to 
the orbit plane. Call this rotated axis set (X", Y", Z"). Then, 

fS’x} r 1 0 0 n/'S'x^ 

S" = S" y = 0 cos(-ORBINC) -sin(-ORBINC) S’y (3.2.1.6) 

S"z J U sin(-ORBINC) cos(-ORBINC) 

Here the sign applied to ORBINC has been chosen to be consistent with the form of the rotation 
matrix in subroutine XROT. 

Now, S" must be rotated by the negative of the orbit position angle p, -RANORB, about 
the Z" axis to align the new X"' axis with the satellite X axis using 

r s "'x\ cos(-RANORB) -sin(-RANORB) (K/S" x \ 

S'" = S’" y = sin(-RANORB) cos(-RANORB) 0 S" y (3.2.1.7) 

\S"'z) 0 0 lJ(S"zJ 

where the sign applied to RANORB has been chosen to be consistent with the form of the rotation 
matrix in subroutine ZROT. 

Last we rename the axes to the corresponding names of the satellite axes: 

X"' -> X 
Y" -> Z 

Z'" Y (3.2. 1.8) 

so that the unit vector to the Sun becomes in satellite coordinates 

S = S’" x i + S"'z j + S’" y k (3.2.1. 9) 

where i, j, and k are unit vectors in the X, Y, and Z directions. Because the distance to the Sun 
from the Earth is large compared to the distance from the center of the Earth to the satellite, parallax 
effects in the coordinate transformations have been ignored. 

Having determined the direction to the Sun in satellite coordinates, it must be determined 
whether the Sun is behind the Earth when viewed from the satellite. It is obvious that the Sun is 
behind the Earth if the distance of closest approach of a ray to the Sun to the center of the Earth is 
less than the radius of the Earth. The task is to determine this distance of closest approach. 
Consider the unit vector S to the Sun (fig. 3.2.1-3). Take the dot product of S and OC, the vector 
from the satellite to the center of the Earth. If this dot product is negative, the Sun is on the 
opposite side of the satellite from the Earth and is unblocked by the Earth. 
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Figure 3.2. 1-3. Sun Ray Intersection With Earth. 


If the dot product is positive, S may intersect the Earth. Noting explicitly that the satellite 
coordinate system is chosen so that die origin is at the satellite and the center of the Earth is at 
(CEARTH, 0, 0), where CEARTH is the negative of the sum of the radius of the Earth RE and the 
satellite altitude, the distance of closest approach is 

TCA = CEARTH sin a (3.2.1.10) 

Noting that 

OC S = CEARTH cos a = CEARTH-S X = RAYCTR (3.2.1.11) 

sin^ a = 1 - cos^ a (3.2.1.12) 

and squaring equation (3.2.1.10) and substituting leads to 

TCA = CEARTH 2 - RAYCTR 2 (3.2.1.13) 

If TCA is greater than RE, the ray to the Sun is unblocked. 
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Ray Intersection With the Earth and Earth-Reflected Solar Radiation. For the 
calculation of Earth-reflected radiation the portion of the Earth visible from the satellite is divided 
into a grid of spherical rectangles of equal azimuth angle and zenith angle extent relative to the 
satellite coordinate system with the X axis as the polar axis. The azimuth grid angles range from 0 
to 360 degrees and the zenith grid angles range from 180 degrees directly under the satellite to 180 
- THETA MX degrees at Earth tangent as seen by the satellite. Here, 

THETAMX = asin (RE/ICEARTHI) (3.2.1.14) 


A ray is generated from the satellite to the angular center of each of the spherical rectangles 
in the grid. If 0 and <)> are the zenith and azimuth angles, respectively, of the angular center of a 

spherical rectangle in the grid, the Cartesian satellite coordinates of a unit vector in the (0,<j>) 
direction is 

RD = cos 0 i + cos <j> sin 0 j + sin <j> sin 0 k = RD X i + RDy j +RD z k (3.2.1.15) 

where, as before, i, j, and k are unit vectors in the X, Y, and Z directions, respectively. 

The ray tracing algorithm given by Glassner (ref.6) has been adapted for determining the 
intersection point of RD with the Earth and the Earth outward normal. Consider the geometry of 
figure 3.2. 1 -4. In this geometry OC = CEARTH i is a ray from the satellite to the center of the 
Earth. The distance of closest approach to the Earth's center along RD is 

TCA = OC-RD = CEARTH RD X (3.2.1.16) 

Because RD is guaranteed to intersect the Earth, TCA is always positive. The half chord distance 
squared is given by 

THC 2 = RE 2 - CEARTH 2 + TCA 2 (3.2.1.17) 

and, therefore, the distance along RD to Earth intersection is 

T = TCA - THC (3.2.1.18) 

Remembering that the satellite is at the origin and that RD is a unit vector, the intersection point 
coordinates are 


PT = RDT = (RD x T,RD y T,RD z T) (3.2.1.19) 

The outward unit normal vector of the Earth at PT is [the center of the Earth is at (CEARTH, 0, 0)] 

NRM = (RD X T - CEARTH)/RE i, RD y T/RE j, RD Z T/RE k (3.2. 1 .20) 

If the dot product of NRM and the vector to the Sun S (equation 3.2. 1.9) is positive; that 
is, 


SUNNRM = S • NRM > 0 (3.2.1.21) 

PT is illuminated by the Sun. The fraction of solar fluence diffusely reflected toward the satellite 
from the spherical rectangle about PT is 
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F = (SUNNRM) (ALBEDO) (dQ/rc) (3.2. 1 .22) 

where SUNNRM is the cosine of the angle between the solar ray and the Earth surface normal at 
the center of the spherical rectangle given by equation (3.2.1.21); ALBEDO is the average Earth 

albedo under the satellite orbital envelope; and d£2 is the solid angle of the spherical rectangle as 
seen by the satellite. 

dO = sin 0 A0 Acp (3.2.1.23) 

where A0 and A4> are, respectively, the zenith angle and azimuth angle extents of the spherical 
rectangle. 

Monte Carlo Scattering. When a ray of solar fluence strikes a surface, it may undergo 
one of three fates: specular reflection, diffuse reflection, or absorption. The probability of these 
fates depends on the surface properties of the surface being struck. The surface property of a fate 
is defined as the probability of that fate occurring. 


A Monte Carlo selection of the fate of the ray is made by choosing a random number 
between 0 and 1. The three surface property probabilities are arranged in the cumulative order 
shown in figure 3.2. 1-5 and the fate is assigned to the bin which contains the random number. 
For example, if the specular reflection is 0.40, diffuse reflection 0.25, absorption 0.35, and the 
random number is 0.53, the ray will be diffusely reflected from the surface. 
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Figure 3.2. 1-5. Cumulative Order of Surface Properties. 

If absorption is selected, the CESH due to the ray is accumulated under the proper heading 
and no further consideration is given the ray. 

If specular reflection is selected, the ray is reflected from the surface as shown in figure 
3.2. 1-6; that is, the component of the ray normal to the surface is reversed. In the figure R is the 
ray to be reflected, N is the surface normal, and S is the specularly reflected ray. R, N, and S are 
all unit vectors. The normal component of R is R*N, so 

S = R - 2(R • N) N (3.2.1.24) 

N 



Figure 3.2. 1-6. Specular Reflection. 

Diffuse reflection is defined as having equal probability for any direction above the surface 
independent of the direction of the incident ray to be reflected. Die reflection direction is randomly 
selected as follows (fig. 3.2. 1-7). In the following all vectors are unit vectors. The surface 
normal vector at the point of ray intersection with the surface is N and T x and Ty are tangent 
vectors. The normal component of the diffusely reflected direction vector D is chosen to be a 
random number r n between 0 and 1 (the normal component of the reflected ray is always parallel to 
the surface normal vector). The components in the T x and Ty directions are chosen as random r x 
and Ty between -1 and +1 . Last, r n , r x , and ry are normalized so that r n ^ + r x ^ +ry^ = l . Then, 

D = r x T x +ry Ty +rn N (3.2.1.25) 
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The reflected ray, either S or D, is ray traced to determine if it hits another surface or not. 
If no surface is struck, the ray is not considered further. If a surface is struck, the ray is checked 
to see if the active side has been struck; if not, the ray is removed from further consideration. If 
the active side has been struck, the point at which the ray strikes the surface is assigned to the 
proper grid point cm the surface and the CESH from the ray is added specular or diffuse, as 
appropriate, to the fluence at that point 

The fate of the ray at this point is Monte Carlo selected and the ray is propagated as 
described above. This process continues until the ray is terminated by absorption, striking the 
inactive side of a surface, or propagation to space, or 100 reflections have taken place. This last 
condition prevents the possibility of an endless loop for rays trapped between perfectly reflecting 
surfaces. 



Figure 3.2. 1-7. Diffuse Reflection. 

Calculation of the solar exposure to a node on the surface requires that the intensity 
normalization of solar rays be discussed. Cumulative equivalent Sun hours (CESH) exposure to a 
flat unshielded surface is the same whether the surface area is large or small or whether one 
considers the exposure to the whole surface or only to a portion of it Therefore, the exposure is a 
fluence; that is, exposure per unit area. A single direct solar ray is given a fluence of SUNTOT / 
ISUNPOS where SUNTOT is the event time in hours and ISUNPOS is the total number of Sun 
positions. A single direct ray of Earth-reflected solar radiation is given a fluence of F (NALB) 
(SUNTOT / ISUNPOS) where F is given by equation (3.2.1.22) and the factor of NALB accounts 
for the fact that Earth-reflected radiation is calculated once every NALB Sun positions. The direct 
exposure to a point on a surface due to a single ray of either direct solar or Earth-reflected solar 
radiation is the fluence of the ray times the magnitude of the cosine of the angle between the surface 
normal at the point and the ray direction. For direct exposure, a single point represents the CESH 
to all points in the nodal area associated with the point. 
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Figure 3.2. 1-8 shows a typical surface divided into six nodes in (£/n) space. In the figure, 

numbers indicate the node number. Direct exposure is calculated at the center of each node in 
space. 



Figure 3.2. 1-8. Nodes and Gridding on a Surface 

A somewhat more sophisticated approach must be taken for calculating scattered CESH 
than for calculating direct CESH. Here, die areas of the nodes on the scattering surface and the 
receiving surface must be taken into account Calculation of the areas of nodes is discussed in the 
Geometric Calculation and Ray Tracing Theory section. Consider figure 3.2. 1-9. In the left half 
of the figure, a large area node scatters a uniform fluence (energy/area, represented by equally 
spaced arrows) toward a small area node. It is clear that only part of the energy scattered from the 
large node is deposited on the small node; however, if the two nodes are parallel to each other, the 
small node receives the same average fluence as the large node reflected. Now, consider the 
opposite case shown in the right half of the figure: a small area node scattering energy to a large 
area node. Here (again, assuming that the two nodes are parallel to each other), all of the energy 
scattered from the small node is deposited on the large node. However, the fluence on the large 
node is smaller than that on the small node because the energy scattered from the small node is 
averaged over the whole area of the large node. The scheme described below properly accounts 
for scattered radiation from one node to another. 



Figure 3.2. 1 -9. Scattering to Different Sized Areas. 
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The energy from direct solar or Earth-reflected radiation falling on a node is the direct 
fluence tim es the projected area, where the projected area is the actual area times the magnitude of 
die cosine of the angle between the ray direction and the node normal. This energy is divided into 
NRAY rays of energy to be scattered from the area. 

NRAY = NINT[(AREA RAY-NRM) MAXRAY / ARMIN] (3.2. 1 .26) 

for direct solar illumination where MAXRAY / ARMIN is the number of rays MAXRAY per 
average nodal area ARMIN and the term in parenthesis is the projected area of the directly 
illuminated node. The NINT function indicates that its argument should be rounded to the nearest 
integer. For Earth-reflected radiation, equation (3.2.1.26) holds if MAXALB is substituted for 
MAXRAY. The energy in each of the NRAY scattered rays is 

SUNTOT ARMIN ^ 7 i 77 ^ 

ESun - ISUNPOS MAXRAY 1 ' * ' ' 


for scattered direct solar exposure and 

F SUNTOT NALB ARMIN 
EEarth - ISUNPOS MAXALB 


(3.2.1.28) 


for scattered Earth-reflected radiation, with F defined in equation (3.2.1.22). 

When one of these rays of scattered energy strikes a surface, the energy in the ray is 
assigned to node which the ray struck. The fluence contributed by this ray is the energy divided by 
the area of the node to which the energy has been assigned. To account for the distribution of 
scattered energy from the node, each ray is scattered from a randomly selected point in the node. 

Geometric Calculation and Ray Tracing Theory. The surface geometry 
description and ray tracing algorithms used in SOLS HAD were developed by Dr. R. C. Corlett 
some years ago. The four geometry and ray tracing subroutines used, RAYG1, RAYG2, 

RAYG3, and RAYG4, are little modified except for style from those originally developed by 
Corlett. The algorithm description given below follows closely the description originally given by 
Corlett (ref. 7). Briefly, the four subroutines perform the following functions. RAYG1 reads the 
parametric description of the surfaces making the object and converts this description to the form 
used internally. RAYG2 returns the Cartesian coordinates of a point specified in internal 
coordinates. RAYG3 returns the normal to a surface at a given point as well as two orthogonal 
vectors tangent to the surface at the point RAYG4 traces a ray in a given direction from a point 
and determines the closest surface which it intersects. 

Surface Geometry Definition. Figures 3.2.1-10 through 3.2.1-14 describe the geometry of 
the five p rimar y surfaces: the trapezoid, cylinder, cone, disk, and sphere, respectively, or portions 
of them. Each surface has an unambiguously defined positive (+) side and negative (-) side and a 

curvilinear coordinate system (£,, T|). The (£, ri) coordinate system is normalized such that the 

primary surface occupies the region [0 < J; £ 1, 0 < r\ < 1 ] exactly. 
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Figure 3.2.1- 1 1. Cylinder Geometry Definition. 
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Figure 3.2.1-14. Sphere Geometry Definition. 

Each primary surface is divided into nodal surfaces. Each nodal surface is bounded by 

constant % and T| boundaries. The positive and negative sides of the nodal surfaces match those of 
their primary surfaces. Figure 3.2.1-15 illustrates the layout of the nodal surfaces for NC = two 

nodes in the ^ direction and NN = three nodes in the T] direction. The primary surfaces are defined 
in physically convenient parameters as described below. 
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Cylinder 

Figure 3.2.1-15. Node Arrangement on Surfaces. 

The trapezoid is described by three points Pi, P 2 , and P 3 oriented as shown in figure 
3 . 2 . 1-10 and the ratio 



(3.2.1.29) 


where 0 < X < 1 and ( p 2 - p .)ll( p 3 ^ 4 ). When Pi, P 2 , and P 3 are oriented in clockwise 

order, the positive side of the trapezoid is visible. This is the normal right hand rule. 

The cone and cylinder sections (figs. 3.2.1-11 and 3.2.1-12) are defined in similar manner 
because the cylinder is a special case of the cone. Point A is the center of the circle of radius R 
(cylinder) or Rmav (cone) at the base of the cylinder or cone and normal to the axis. Vector E is 
parallel to and co-directed with the axis of the cone, h is the ratio the height H to R (cylinder) or 
Rmax (cone). Reference point P re f defines azimuthal angle 0 = 0 and 0 is measured in the right 
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hand direction about the axis vector E. 0m is the maximum value of 0. For the cone R min = g 
Rmav where g is a dimensionless constant between 0 and 1. From this we see that the cylinder is a 
cone with g = 1. 

The disk (fig. 3.2.1-13) is defined similarly to the cone with h deleted, Rmin = ri and 6 
measured in the left hand direction about the axis vector E. 

The sphere (fig. 3.2.1-14) is defined in similar manner to the cone with g and h deleted and 
bounding polar angles <|>i and <J>2 added to define the polar extent of the sphere section. Vector E is 
co-dirccted with the half line <J> = 0. 

A fundamental reference coordinate system for each primary surface is defined as follows. 
For the trapezoid. 


( P 2- P l)*( P 3- P l) 

‘-|( p 2 - p ,)*( p 3 - p ,)| 


*1 


( P 2 ~ P l) 

K P 2- P l)| 


y„ = *xy, 


and for the cylinder, cone, disk, and sphere. 


(3.2.1.30) 


(3.2.1.31) 

(3.2.1.32) 



EX ( P ref- A ) 

y2 "| lt>< ( P ,. f - A ) 


(3.2.1.33) 


(3.2.1.34) 


Y 1 =y 2 X 


(3.2.1.35) 


To avoid numerical computation problems, the following restrictions are placed on the 
above definitions. For the trapezoid, the included angles between sides must be greater than ca, 
where ca is set to 0.1 degree in RAYG1. This requirement prevents two adjacent sides of the 
trapezoid freon being parallel. 


For the other primary surfaces, the magnitude of the angle included by (Pref - A) and E 
must be less than ca, where ca is as defined above. Other variables used to define cylinders, 
cones, disks, and spheres are restricted to the following ranges: 
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O<0 m < 360 deg, 

0<g< 1, 

0 < < <t»2 < 180 deg, 

<t>2 - 4>1 > o. 

Surface points on each type of primary surface arc defined in terms of input parameters and 
the dimensionless coordinates ^andr|. For the trapezoid a point P is located at 

P= P, + n( p 2 - P 0 + 5[P 4 + H( P 3 - P 4 ) - P , ( P 2 - P 0] 


= p , + S( p 4 - p,) + n ( p 2 - p ,) ( p , - p 2 + p 3 ' p 4 ) 

For the cylinder, cone, and disk 

(3.2.1.36) 

P= A + ^rh* + r y [1+^(g- 1)] 

where g — 1 for a cylinder and h = 0 for a disk. For the sphere. 

(3.2.1.37) 

P = A + r [*cos<|> + y si n<J>] 
In the above two equations. 

(3.2.1.38) 

0 = T l e max 

(3.2.1.39) 


(3.2.1.40) 

y = y. | COS0+ y g Sin6 

(3.2.1.41) 

Differential areas cL4 for the five primary surface types may be expressed in terms of input 
parameters and E, and tj as follows. For the trapezoid 

dA = [(P 2 - Pi) * yil [(P3 - Pi) ' y 2 l [1 + ft -1) £1 d£ dq 
For the cylinder, cone, and disk 

(3.2.1.42) 

dA = R 2 0m [h 2 + (g - l) 2 ] 172 [1 + (g-D ^ d£ dr) 

with g = 1 for the cylinder and h = 0 for the disk. For the sphere 

(3.2.1.43) 

dA = R 2 0 m (<t>2 ~ <t>l) sin [<j>l + £(<l>2 " $1)1 d£ dq 

(3.2.1.44) 


Because an individual nodal surface is bounded by < % < %2 and H 1 < “H < T 12. the areas 
of nodal surfaces can be determined by integrating 
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(3.2.1.45) 


S2 TJ2 

A = JdS JdqdA(S.q) 

Si ill 

which yields for the trapezoid 

A = [(P2 - Pi) • yil [(P3 - Pi) * y2l (T12 - Til) (S2 ~ Si) [1 + (S2 + Si) ft -l)/2] (3.2.1.46) 

For the cylinder, cone, and disk 

A = R2 0 m [h2 +( g -i)2]l/2 (T12 . q t ) (^ 2 - Si) [1 + (S2 + Sl) (g -l)/2] (3.2.1.47) 

with, as before, g = 1 for the cylinder and h = 0 for the disk. For the sphere 

A = R 2 e m (q 2 -qi) {cos [<j>i + Si (<t>2-<t>l)]-cos[<|)i+S2(<l>2-<t>l)]} (3.2.1.48) 

Note that in equations (3.2. 1.46) through (3.2.1.48) that the areas are direcdy proportional to q, 
but are nonlinear in S- Further, the expressions for area may be separated into three factors: a 
constant factor, the factor (r|2 - tu), and a S dependent factor f(S2» Si)- Thus, if the area of a 
node is known, the area of a portion of the node Ap bounded by H 1 < T| <Tip, qp< 712 , and Si < 
S < S2 is given by 


Ap=A (Tip -TU)/ 012 -m) (3.2.1.49) 

Similarly if 4 p bounded by q i < q < q 2 and Si < S < Sp» Sp < S2. 

A V =A f(Sp, Si) / f(S2, Si) (3.2.1.50) 


A nodal surface is bounded by Si < S < S2 and q l < q < q 2 . It is desirable to determine 
the (S, T|) position of a point given its relative position Rj and Rk on a nodal surface, where 0 < Rj 
< 1 is in the q direction and 0 < Rk ^ 1 is in the S direction. 


R J = ( t i- T ti)/( T i 2 -^,) or n-n, + R j (’i 2 -n 1 ) 

and 

[i+(p-i>a 2 -[i+(P-i)?J 

R k — — "T 

[i + (p-i)S 2 j -[ 1 + (P-i)5J 


or 


(3.2.1.51) 


(3.2.1.52) 
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(i- 7 d * (P ■ - -0 jj + R ,( t, -tQ(P- Dp + <P- 1 >( t g -tfl) 

1- p 


i f p= 1 


i f p * 1 

( 3 . 2 . 1 . 53 ) 


For the trapezoid, p = X; for the cylinder, p = 1; for the cone and disk, p = g. 


For the sphere, 

cosfy + 5, (* 2 - ♦,)]- cosfy + § (<t> 2 - ♦,)] 
‘ k= co^o, + 4, (♦, - ♦,)]- co{*, + ♦,)] 


( 3 . 2 . 1 . 54 ) 


or 

% = [c0S" 1 (C0S[<|> 1 + ^(*2 - < > 1 )] + 

+ ^2 (*2 " * 1 )] " C0! {*2 + ^(*2 ~ * l )]}) 

_<t> J (* 2 ”* l ) ( 3 . 2 . 1 . 55 ) 

To avoid the inverse cosine operation in the above equation and to simplify other calculations, it is 
convenient to rewrite the equation using the following change of variable 


cos(j> 1 - cos[<j> 1 + $(<|> 2 - 4^)] 
cos(<|> 1 ) — cos(<(> 2 ) 




Thai 

cos<|> = cos<|> 1 +£(cos(<|> 2 ) - cos^)) 


( 3 . 2 . 1 . 56 ) 


( 3 . 2 . 1 . 57 ) 


which leads to the following simplified equation for the sphere section which will be used hereafter 


^“^1 + R k (^2 *» i ) ( 3 . 2 . 1 . 58 ) 

Generation of Local Unit Normal and Tangent Vectors. The unit normal vector N directed 

outward from the positive side of a nodal surface at point P corresponding to (!;, rj) on some 
primary surface and the orthogonal unit tangent vectors T i and T2 are defined such that 

T 1 xt 2 = N (3.2.1.59) 

Then, for the trapezoid and disk 
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N=z, 


T i = y r 


T 2 = Y 2 


(3.2.1.60) 


For the cone and cylinder, define 


cosa 


sina 


= h/ 1 /h 1 2 + (g-1) i 

- 0-9)1 J' 


h 2 + (g-1) 2 


(3.2.1.61) 


(3.2.1.62) 


where g = 1 for the cylinder and g is the ratio of the minimum to the maximum radius for the cone. 
Then, for the cone and the cylinder 


N = y cosa + x sina 


t.| = — y sina + z cosa 


T 2 = KXT ( = — y 2 COS 0 + y 1 sine 

and for the sphere 
N= z COS<|> +y sin<|> 

T 2 = -y 2 COS0 + y 1 Sine 
Ti=T 2 xN = .Sin4 > -y COS4» 


(3.2.1.63) 

(3.2.1.64) 

(3.2.1.65) 

(3.2.1.66) 

(3.2.1.67) 

(3.2.1.68) 


In the above equations 6, <(>, and y are defined in equations (3.2.1.39) through (3.2.1.41). 

Rav Tracing. This section describes the method for determining whether a ray in a given 
direction from a point on a nodal surface intersects another surface, and, if so, whether that surface 
is the nearest surface intersected. Suppose Lm is the shortest ray length previously found from the 

point of origin P 0 and that the ray is co-directed with unit vector v to an intercept on some nodal 
surface. The problem is to consider some new primary surface and answer the following 
questions: 

1. Does this ray intercept the primitive surface 6 containing primary surface of concern with 
positive intercept distance L < Lm? 

2. If so, is the intercept point within the primary surface? 

3. If so, where within a nodal surface is the intercept point contained? 


^The primitive surface of a trapezoid or a disk section is the infinite plane which contains it. The primitive surface 

of a cylinder section is the complete infinite length cylinder. The primitive surface of a cone section is the complete 

infinite cone. The primitive surface of a sphere section is the complete sphere. 


74 



If the primitive surface is quadratic, that is, a cylinder, cone, or sphere, there will in general be two 
intercept points.The nearer one (the one with the shorter L) is considered first until one of the 
above questions has a negative answer, after which the sequence of questions is asked of the 
further point If the first two questions have positive answers, die point of intercept is taken as 
interim ray termination point P and Lm is replaced by L and the nodal surface within the primary 
surface is calculated. If the first question has a negative answer, the remaining questions need not 
be considered. The process is repeated until all primary surfaces have been considered. 

The three questions are considered separately in the following subsections. 


Does the rav intercept the primitive surface with 0 < L < Lr rj? The intercept point P may be 
expressed 


P= P 0 + Lv 


(3.2.1.69) 


Using equations (3.2.1.36) through (3.2.1.38) as appropriate to express P, this vector equation 
contains three scalar unknowns, L, 2;, and rj. In dealing with this equation, it is convenient to 
introduce the quantity B, 

B=Pi-P 0 (3.2.1.70) 

for die trapezoid or 

B = A - P 0 (3.2.1.71) 


for cylinder, cone, disk, or sphere sections. 

For the trapezoid and disk section primary surface types, taking the dot product of each 
side of equation (3.2. 1.69) with the vector z gives directly 

L(v x) = (B • x) (3.2.1.72) 


from which sign (B*z) = sign (vz) and the criterion 0 < L < Lm may be written as 

0 < IB*zl < Lm Iv-zl (3.2.1.73) 

For the cylinder section or cone section primary surface types the calculations are more 
complex. Taking the dot product of each side of equation (3.2.1.69) with the vector z yields 

^ r h =L (v*z) - (B*z) (3.2.1.74) 

Then, after rearranging equation (3.2.1.69) so that the term r y [1 + £ (g-1)] is alone on one side, 
squaring each side, and substituting for ^ using equation (3.2.1.74), there results 


AL 2 -2BL-C = 0 


(3.2.1.75) 


where 
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(3.2.1.76) 


A = 1 - (v-e) 2 


B = (V’B) - (v*e) [(B*e) - 5] 

(3.2.1.77) 

C = (B-e) [(B-e) - 2 5] + r 2 - B-B 

(3.2.1.78) 

e = z(l+ 7 2 )!/ 2 

(3.2.1.79) 

8 = ry(l+y 2 )!/ 2 

(3.2.1.80) 

Y= (g - l)/h 

(3.2.1.81) 


and r = R for the cylinder and Rmax for the cone. 

For the sphere section primary surface type, after rearranging equation (3.2. 1.69) so that 
the term R [z cos <]> + y sin <|)] is alone on one side and squaring each side, there results similarly 

AL 2 -2BL-C = 0 (3.2.1.82) 

where A = 1, B = (v*B), C = r 2 - B*B. 


Thus, for all three quadratic primary surface types 
L A = B ± VD 


(3.2.1.83) 


with D = B 2 + A C. 

Note that the vector e and the scalar 8 are independent of ray direction or point or origin. In 
stating a criterion for the existence of L such that 0 < L < L®, it is convenient to take A positive. 

The case A = 0 may be disregarded because its probability is negligibly small with v chosen 
randomly. If A < 0, the quadratic equation (3.2.1.83) for L is preserved if the sign of each of the 
three coefficients A, B, and C is arbitrarily reversed, which will be assumed to have been done if 
necessary to make A positive. 

Is P within the primary surface? In this subsection it is assumed that there does exist L 

such that 0 < L < Lm and that L has been calculated from equation (3.2.1.72), (3.2.1.74), or 
(3.2.1.83) as appropriate for the primary surface type. For all primary surface types but the disk 
section, the first step is to calculate the normalized coordinate t, directly and to test for 0 < £ <1. 
Only for the disk section primary surface type do computational savings result from an indirect 
calculation. 

For the trapezoid primary surface type, taking the dot product of each side of equation 
(3.2.1.69) with unit vector y2 gives 

£ = w • (L v - B) (3.2.1.84) 
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where w = y 2 / [y 2 ‘ (P4 - Pi)]- 

For the cylinder and cone section primary surface types, equation (3.2. 1.74) can be 
rewritten as 

£ = w*(Lv-B) (3.2.1.85) 

where w = z / (r h). 

For the disk section primary surface type, inspection of (3.2.1.69) (noting that h = 0) 
shows that 

ILv-B! = R[l +£(g-l)] (3.2.1.86) 

or 

5 = [l-ILv-BI/R]/(l-g) (3.2.1.87) 

For this primary surface type the criterion 0 < £ < 1 is most conveniently applied in the form 

R2 > L 2 + B • B - 2 L (v - B) > R 2 g 2 . (3.2.1.88) 

For the sphere section primary surface type, taking the dot product of each side of equation 
(3.2.1.69) with unit vector z yields 

cos <j> = z*(Lv-B)/R (3.2.1.89) 

which, upon use of equation (3.2.1.57), becomes 

£=v-w*(Lv-B) (3.2.1.90) 

where w = z /[R (cos 4»i -cos <j>2)] and v = cos <(>i / (cos 4>1 -cos <|> 2 ). 

If 0 < £ < 1, then the second step, determining if 0 < 1 ) < 1, is carried out. It is assumed 
that both L and £ have been evaluated. 

For the trapezoid primary surface type, taking the dot product of each side of equation 
( 3 . 2 . 1 . 57 ) with yi and noting the definition of the length ratio X yields 

ti[1+£(X-1)] = u*(Lv-B)-v£ (3.2.1.91) 

where u = yi / [yi • (P 2 - Pi)] and v = yi • (P 4 - Pi) / [y 1 - (?2 - Pl)l- 

For the other four primary surface types, equation (3.2.1.69) is solved for cos 0 through 
dot multiplication by yi and for sin 0 through dot multiplication by y 2 - The solutions are 

[1 + S(g- l)]cos0 = Yl-(Lv-B) (3.2.1.92) 
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[1 + £ (g - 1)] sin 0 = Y 2 • (L v - B) 


(3.2.1.93) 


where Yi = y l / r and Y 2 = y2 / r and r = R or R m ax depending on surface type. Then, 0 < r\ < 1 
if sin 0m ^ 0 and sin 0 > 0 and cos 0 > cos 0m or if sin 0 m < 0 and sin 0 > 0 or cos 0 < cos 0 m . 

Within what nodal surface is P? It is assumed in this subsection that \ and t| have already 
been calculated and that the number of nodes in the £ direction NC and the number of nodes in the 
T| direction NN for the primary surface are known. The primary surface is a square in (£, r|) 
space, that is, 0 < ^ < 1 and 0 < q < 1, and each nodal surface is a rectangle in (%, q) space 

bounded by 0 ^ %\ < £ < £ 2 < 1 and 0<T)l<T|<T| 2 £l. Each nodal surface on a primary 
surface is given an index number. Figure 3.2.1-15 shows the arrangement of nodal surfaces on 

each of the five primary surfaces for NC = 2 and NN = 3. The position of £ in terms of nodal 
surfaces within the primary surface is given by 

N£ = (£ NC) (3.2.1.94) 

where the term in parenthesis is rounded up to the next integer. The position of q in terms of nodal 
surfaces within the primary surface is given by 

Nq = (q NN) (3.2.1.95) 

where again the term in parenthesis is rounded up to the next integer. Now, given the nodal 
surface coordinates as N£ and Nq, the index IN of the nodal surface containing £ and q is given 
by 

IN = Nq + NN (N£ - 1). (3.2.1.96) 

Two-Dimensional Interpolation. CESH are calculated at the (£,q) centers of the 
nodes on surfaces. However, for proper display during plotting, CESH are needed near the 

comers of the nodes 7 . The following scheme interpolates CESH from the (^.q) centers of the 

nodes to the comers of the nodes. It is most convenient to perform the interpolation in (£,q) space 
rather than in Cartesian (x,y,z) coordinates. 

As noted previously, each surface is a square in (£,q) space with 0 < 2; < 1 and 0 < rj < 1. 
For a surface divided into NC equal nodes in the ^ direction and NN equal nodes in the rj 
direction, the dimensions of each nodes are 1/NC in the £ direction and 1/NN in the T] direction. 
The node centers are at ( c Ej, °r|i) = ([j-0.5]/NC, [i-0.5]/NN) with j = 1,2,...,NC and i = 

1,2,...,NN. Figure 3.2.1-16 shows a typical surface representation in (£,q) space with NC = 2 
and NN = 3. The figure also shows the location of the points near the comers of the nodes, whose 

?To avoid ambiguity in node assignment, the points "near the comers" are chosen to be set in from 
the comers by 1% of the node dimension. Hereafter, references to the comers of nodes will be 
synonymous with "near the comers." 
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values of CESH are to be interpolated. Inspection of the figure shows that only two of the points 
whose values are to be interpolated lie within the perimeter of the points at the centers of the nodes. 
Because extrapolation is, in general, less accurate than interpolation, die edges of open surfaces 8 
are padded with extra rows and columns as follows: The lower edge is padded by repeating the 

CESH values for c £j below the surface and the top edge by repeating the CESH values for C ^NC 

above the surface. The left and right edges are padded using the CESH values for °ni and °nNN, 
respectively; and the comers by repeating the CESH values at the corresponding node centers at the 
comers of the surface. Closed surfaces are padded in similar manner except that the left column is 

padded with CESH values of the °T|NN column and the right column with CESH values from the 
Ct|i column. The centers of the nodes added for padding are at ( c E,j, °n = -0.5/NN) and ( c ^j, °n = 
1 + 0.5/NN) for the left and right columns, respectively, and ( c £ = -0.5/NC, °r|i) and ( c £, =1 + 
0.5/NC, °r|i) for the bottom and top rows, respectively. 

These preliminaries concluded, the description of the interpolation algorithm may begin. 
Suppose one wishes to interpolate the the value of CESH y(^,Tl) at (^,1)), given that c ^j < ^ ^ 
c £,j+l and °ru < r\ < Cqi+i. The interpolation is performed in two steps. First, one-dimensional 
interpolations are performed on the rows for each row, 1 < j ^ NC, to give the one-dimensional 
array of CESH ( c £jlr|) where ( c ^jlq) indicates that Ti is fixed. Then, this array is interpolated to 
give the value of CESH at 

The interpolation scheme used is Stineman's consistently well-behaved method of 
interpolation (ref.4). The following explanation of the method quotes parts of reference 4 directly 
except for changes of equation and figure numbers to be consistent with this document and 
additions and deletions specific to the problem addressed by this document. 

The interpolation scheme has the following three properties. First, if values of the 

ordinates of adjacent specified points (the CESH at the node centers in the ^ ort) direction) change 
monotonically and the slopes of the line segments joining the points change monotonically, then 
the interpolating curve and its slope will change monotonically. Second, if the slopes of the tine 
segments joining the adjacent points change monotonically, then the slope of the interpolating 
curve will change monotonically. Third, suppose that the first two conditions are satisfied by a set 
of points, but that a small change in the ordinate or slope at one of the points will result in one of 
the conditions no longer being satisfied. Then, making this small change in the ordinate or slope at 
a point will cause no more than a small change in the interpolating curve. 

The interpolation scheme is implemented as follows. Assume that the ordered sequence of 
data (xj, yj), j = 1, 2, ..., n is given where (xj, yj) are the rectangular coordinates of the j^ 1 point 

on the curve (here xj may be either c E,j or °tu and n may be NC + 2 or NN + 2 as appropriate and 
yj is the corresponding CESH). Let yj’ be the slope of the curve at the j* point and require that xj 
< xj+i for j = 1, 2, ..., n-1. Calculation of yj' will be considered presently. Before doing so, it is 
noted that the yj' are calculated most accurately if the xj and yj have roughly equal ranges. To 


8 An open surface is defined as one such that the Cartesian coordinates of (^,p = 0) and (£,r| = 1) 
are different, for example, a trapezoid. A closed surface is defined as one such that the Cartesian 

coordinates of (£,T| = 0) and (£,1) = 1) are identical; for example, an annulus of a cylinder or of a 
sphere. 
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ensure this, the values of c £,j, °rii and CESH( c ^j, Crjj) are each scaled to be in the range 0 to 1 
before interpolation begins. After completion of interpolation, the interpolated y are rescaled back 
to their true values. 



T) 


+ Points at centers of nodes 
• Points near corners of nodes 

— Node boundaries on surface 

— — Rows and columns of nodes padded onto surface perimeter 

Figure 3.2.1-16. Surface in Space With NC = 2, NN = 3. 


Given x such that xj < x < xj+i, the procedure for calculating y (the corresponding 
interpolated value) is the following. The slope of the line segment joining the two points is 
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(3.2.1.97) 


Sl . = m i - y j 

J Xj+1-Xj 

On the line segment the ordinate corresponding to x is 
yo = yj + sj (x- xj) (3.2. 1.98) 

Next, 

Ayj = yj + yj’ (x - xj) - y 0 (3-2. 1 .99) 

where Ayj is the vertical distance from the point (x, y 0 ) to a line through (xj, yj) with slope yj', as 
shown in figure 3.2.1-17. Similarly, 

Ayj+l = yj+l + yj+l' (x-xj+i)-y 0 (3.2.1.100) 

is the vertical distance from the point (x, yo) to a line through (xj+i, yj+l) with slope yj+l', also 
shown in figure 3.2.1-17. The product Ayj Ayj+l is then calculated and tested. 



Xj X x j+1 

Figure 3.2.1-17. Consistently Well-Behaved Interpolation for Ayj Ayj+l > 0. 
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If yj’ = sj, then the line through point (xj, yj) will coincide with the line segment joining 
points (xj, yj) and (xj+i, yj+i), and Ayj = 0. Similarly, if yj+i' = sj, then Ayj+i = 0. If either or 
both Ayj or Ayj+i are zero, then the product Ayj Ayj+i = 0 and y is simply 

y=yo (3.2.1.101) 


If Ayj Ayj+i > 0, then (as in fig. 3.2.1-17) Ayj and Ayj+i have the same sign, and 


y = yo 


+ _Ayj_Ayj ± i_ 
Ayj + Ayj+i 


(3.2.1.102) 


Equation (3.2.1.102) always determines the point (x, y) inside the triangle UK of figure 3.2.1-17. 
The slope of the interpolating curve matches the given slopes at the given points. The slope 
changes monotonically between xj and xj+i, so the interpolating curve is always concave toward 
the line segment joining the two points. 


If Ayj Ayj+i < 0, then the geometry is as in figure 3.2.1-18 and there must be an inflection 
point between xj and xj+j. In this case, 


y = y , -AyiAyj+iJx-- xj±x - x,+i) 

(Ayj - Ayj+i) (xj+i - xj) 


(3.2.1.103) 


Equation (3.2.1.103) always determines the point (x, y) inside the quadrilateral JIKL of figure 
3.2.1-18, where the vertical distance LO equals the vertical distance OI. The slope of the 
interpolating curve matches the slopes at the given points xj and xj+i. The interpolating curve 
intersects the line segment JK at its midpoint 


The rationale for equation (3.2.1.103) may be understood by considering the case where yj 
is significantly greater than sj, the slope of line segment JK, but yj+i is nearly equal to sj (fig. 
3.2.1-17 or 3.2. 1-18). Regardless of whether yj+i is greater or less than sj, points I and L will 
be very close to point J and the interpolating curve will be very close to line segment JK. Thus, a 
change of yj+i from slightly more than sj to slightly less than sj will cause only a slight change in 
the interpolating curve. This illustrates the third property, described above, of the interpolation 
scheme . 


The derivatives yj’ are calculated by the following procedure, which ensures that they have 
the properties required by the interpolation method. In figure 3.2.1-19 let I, J, and K be any three 
consecutive points. Point J may be above or below the line segment joining I and K as shown in 
figures 3.2.1-19a and 3.2.1-19b, respectively. The requirements on the derivatives are satisfied if 
yj' has a value between the slopes of the line segments U and JK. That is, for figure 3.2.1-19a, it 
is necessary that slope(U) > yj' > slope(JK), while for figure 3.2.1- 19b, slope(U) < yj* < 
slope(JK). Another point is that if, for example, line segment U is much shorter then JK, it may 
easily be seen that a smoother overall interpolation curve will result if yj' is nearly equal to the 
slope of U. 
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Figure 3.2.1-18. Consistently Well-Behaved Interpolation for Ayj Ayj+i < 0. 


All of the above considerations are satisfied by setting yj’ equal to the slope at point J of a 
circle through points I and K. Thus, 

(yj - yj) ((xk - xj ) 2 1 (yk - yi) 2 ) ± (yk - yj) ((*j - *i ) 2 ± i y\ - (3.2.1.104) 

yj “(xj - xi) ((xk - xj) 2 + (yk - yj) 2 ) + (xk - xj) ((xj - x0 2 + (yj - yi) 2 ) 

Equation (3.2.1.104) takes care of all interior points. The end points require special attention. 
From the requirements on the yj' it may be seen that the slope at an end point must have the same 
sign as the line segment from the end point to the next point In line with the use of equation 
(3.2.1.104), one might draw a circle through the first (or last) three points. However, if either 
point I or K of figure 3.2. 1- 19 is an end point, the particular arrangement of the points causes the 
slope of the circle at I or K to have the wrong sign. The sign of the slope is not necessarily wrong 
at the end point, but it may be wrong. 
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The problem divides into two cases. To simplify the notation, let M and subscript m 
designate either point I or K, whichever is an end point, and let s be the slope of the line segment 
joining points J and M. The first case occurs when s is steeper than yj'. In this case, a parabola 
through J and M with slope yj' at J has a slope at M which meets the requirements and is 
reasonable. Stated more precisely, if s > 0 and s > yj' or if s < 0 and s < yj', then 

ym' = s + (s- yj’) = 2s - yj' (3.2.1.105) 
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The second case occurs if neither condition for the first case is satisfied. In this case, the 
term in parenthesis in equation (3.2.1.105) is multiplied by a factor between zero and one which 
assures that ym' is always the same sign as s. The result is 


Isl (s - yj') 
ym s + Isl + Is - yi'l 


(3.2.1.106) 


3.2.2 MDDB Version 2.0 Algorithm 

The primary mode of operation of MDDB is one of interactive prompt and user response to 
the prompts. Major computational routines to specify the positions of points on surfaces for 
TECPLOT display are CONNECM, RAYG1M, RAYG2M, and RAYG3M. These are modified 
versions of subroutines of similar name in SHADOW (ref. 1) and in SOLSHAD. Users interested 
in the algorithms used in these routines should refer to the algorithm descriptions in section 3.2.1. 
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APPENDIX A: SAMPLE PROGRAM RUN 


Appendix A contains a series of sample files that have been constructed by taking a sample 
geometry and running it through the entire process to its final output The sample geometry 
structure was devised so that each of the five primary surface types is used once and so that, when 
SOLSHAD is run, the effects of multiple reflections are apparent The active side of the sphere, 
cone, and cylinder are on the outside. The active side of the disk and the trapezoid are toward the 
center of the structure. The surfaces have a variety of different surface properties. The ram 
direction (direction of the satellite velocity vector) is along the positive z axis; that is, from the 
center of the disk and pointing through the point of the cone. The orbit is chosen to be a 28.5 deg. 
inclination and a 450 km altitude, with all possible random sun and satellite positions. 

A.1 3D Graph of Surfaces 

Figure A1 is a black and white image of an MDDB TECPLOT color display of the sample 
file geometry. (The structure has been rotated from its original position so that it can be viewed 
more clearly.) An MDDB TECPLOT display will be colored as labeled in figure A1 . The active 
side of each surface is indicated by a surface normal direction arrow. The direction arrow for the 
blue square is not visible because it points down into the sphere. 

When the structure is viewed from the top, the intersection of the sphere and square is an 
example of where TECPLOT sometimes draws lines where it is not supposed to. (See section 
2.2.3 about problems with TECPLOT.) However, this is not shown in Figure Al. 


A1 



Figure Al. Black and White Image of the Sample 
Geometry as Viewed Using MDDB and TECPLOT. 


A.2 Input File Generated by MDDB 

This is the file generated by MDDB. Note that the cone ratio is taken to be 0.001 instead of 
0.000. (See sec. 2.2.3 about problems with TECPLOT.) 

File on PC: sample.sol 

File cm convex: solmicro/solshad.in_sample 

SOLSHAD sample input with 5 surfaces 
1 5 5 T T T 

-4.0000 -4.0000 

-4.0000 4.0000 

4.0000 4.0000 

PLANE 


15.1000 

15.1000 

15.1000 

1.0000 
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PLANE MATERIAL 



.5000 .5000 

.0000 

SPEC, DIFFUSE, ABSORPTION 

2 5 24 T T T 

.0000 


.0000 1.0000 

.0000 


.0000 .0000 

1.0000 


.0000 .0000 

3.0000 


1.6667 

CYLINDER 

CYLINDER MATERIAL 

.2500 .7000 

.0500 

SPEC, DIFFUSE, ABSORPTION 

3 6 24 T T T 

.0000 


.0000 1.0000 

.0000 


.0000 5.0000 

1.0000 


.0000 .0000 

3.0000 


1.6667 .0010 

CONE 

CONE MATERIAL 

.3000 .6000 

.1000 

SPEC, DIFFUSE, ABSORPTION 

4 4 24 T T T 

.0000 


.0000 1.0000 

.0000 


.0000 .0000 

1.0000 


.0000 .0000 

6.0000 


.5000 

DISK 

DISK MATERIAL 

.0500 ,9000 

.0500 

SPEC, DIFFUSE, ABSORPTION 

5 6 24 T T T 

.0000 


.0000 1.0000 

.0000 


.0000 12.5000 

1.0000 


.0000 .0000 

2.5000 


.0000 180.0000 

SPHERE 

SPHERE MATERIAL 

.0500 .9000 

.0500 

SPEC, DIFFUSE, ABSORPTION 

0 0 0 F F F 

END OF GEOMETRY 

1994 6 19 0 

.0000 

Y M D H START 

1994 10 27 12 

.0000 

Y M D H END 

.0000 360.0000 


ORBPOS (1) , ORPBOS (2 ) 

.0000 360.0000 

F 

PRECES ( 1 ) , PRECES ( 2 ) , TIMCOR 

28.5000 


ORBINC 

.2460 6 . 3780E+03 

4 . 5000E+02 

ALBEDO, RE, HS AT 

10 5 10 

1000 T F F MAX RAY, MAX A LB, NALB, ISUNPOS, SHORTL, CHECK 

5 5 


NTHETA, NPHI 

997531 

ISEED 



360.0000 


360.0000 


360.0000 


360.0000 


INACT 
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This is the file that is generated by SOLSHAD as standard output, which is generally redirected to a file. 
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AREA OF AVERAGE NODE - 7.1185E-01 

THIS NODE MAY HAVE UP TO 10 RAYS DIRECT SOLARS C ATT ER ED FROM IT 

AND UP TO 5 EARTH REFLECTED RAYS SCATTERED FROM IT 

PER POINT ON EARTH 
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TITLE-" SOLSHAD aample Input with 5 surfaces" 

VARIABLES-NODE, PRIM, SPEC, DIFF, ABSORPTION, TOTAL, PRIMA, SPECA, DI FFA, ABSA, TOTAL, GRTOTAL, X, Y, 2 
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Figure A3. An X- Y Graph of the Solar Exposure as a Function of Distance Along 
the Diagonal Line of the Square Shown in Figure A2. 




APPENDIX B: 


FORMAT OF OUTPUT FILES FOR TECPLOT 


B.l TECPLOT Output File From SOLSHAD 

SOLSHAD produces an ASCII file of the nodes, solar exposures, and connectivity matrix 
formatted for processing by PREPLOT to create a binary input file for TECPLOT. TECPLOT 
displays this file as a plot of the object with the fluxes color coded on the surfaces. This file is 
namwi TAPE7 and is opened in SOLSHAD (subroutine CONNEC) as STATUS = 

'UNKNOWN', FORM = FORMATTED’. STATUS = UNKNOWN’ causes any existing 
version of TAPE7 to be overwritten. The last line of the batch file run.shadow gives the file 
TAPE7 a new name. Therefore, users should take care to change the name of the file in 
run.shadow if they wish to preserve results of previous SOLSHAD runs. 


Record El Format (lX, , Tn’LE=" , ,A, ,,tl ) 


Entry 


Type 

Variable 

Description of geometry from record B 1 

9-69 

character 

HEADER 


Record E2 Format (' VARIABLES=NODE, PRIM, SPEC, DIFF, ABSORPTION, 
TOTAL, PRIMA, SPEC A,DIFFA,ABS A, TOT AL,GRTOTAL,X,Y,Z') 

This record provides variable names for the columns of data in record E4. 


Record E3 Format (' ZONE T="SURF 1' 

M=M4, , J= , ,I4,'F=FEP( 

DINT') 


Columns 

Type 

Variable 

Number of nodes (points) on the surfaces of the 
structure. 


integer 

NODETOT 

Number of elements in the connectivity matrix. 



NCON 


Record E4 Format (I6,11F10.2,3F10.3) 

This record is repeated NODETOT times, once for the data at each node. 


Entry 

Columns 


Node number. 

1-6 


■■9HI 

Direct solar exposure (CESH) at node. 

7-16 

real 

BaSfflggi 

Solar exposure due to specular reflection of direct 
solar radiation (CESH) at node. 

mim 

WBM 


Solar exposure due to diffuse reflection of direct 
solar radiation (CESH) at node. 

1HKI 

real 


Solar exposure due to absorption of direct solar 
radiation (CESH) at node. 

37-46 

mBtm 


Total solar exposure due to direct solar exposure + 
specular reflection of direct solar radiation + 
diffuse reflection of direct solar radiation (CESH) 
at node. 

47-56 

real 

FLUX(5JL) 

1 EE5SB 55 5* ESHSE5 55 * 1 32L i lEEsSiaMM 

BKIEisMi 

real 


Solar exposure due to specular reflection of Earth- 
reflected solar radiation (CESH) at node. 

1 

real 


Solar exposure due to diffuse reflection of Earth- 
reflected solar radiation (CESH) at node. 


mm 
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Solar exposure due to absorption of Earth-reflected 
solar radiation (CESH) at node. 

87-96 

real 

FLUXA(4,L) 

Total solar exposure due to Earth-reflected solar 
exposure + specular reflection of Earth-reflected 
solar radiation + diffuse reflection of Earth- 
reflected solar radiation (CESH) at node. 

97-106 

real 

FLUXA(5,L) 

Grand total solar exposure (CESH) at node = 
FLUX(5,L) +FLUXA(5,L) 

107-116 

real 

FLUXT(L) 





X coordinate of node 

117-126 

real 

C$AV(1,L) I 

Y coordinate of node 

127-136 

real 

CSAV(2U) 

Z coordinate of node 


real 

CSAV(3,L) 


Record E5 Format (4110) 

This record is repeated NCON times, once for each element of the connectivity matrix. 


Entry 

Columns 

Type 

Variable 

The indices of the four nodes defining an element 
of the connectivity matrix. When the four nodes 
are connected in order, they form a quadrilateral (or 
a triangle if two adjacent points are identical). 

1-10, 11-20, 
21-30, 31-40 


ICONEC 


B.2 TECPLOT Output Files from MDDB 

TECPLOT Binary Output File From MDDB. MDDB produces a binary TECPLOT 
file whenever a TECPLOT display of the object is selected. This file is readable directly by the PC 
version of TECPLOT 5.0x without need for preprocessing by PREPLOT. The binary file is 
named MDDBJPLT and is opened in subroutine GENSURF with a FORTRAN open statement as 
unit 9, STATUS = ’UNKNOWN, FORM = ’BINARY’. STATUS = UNKNOWN causes any 
existing version of MDDB.PLT to be overwritten. The new version of the file is written in 
subroutine CONNECM. Because the file is binary, its format is described by contents of 
individual words rather than by record content. Each word is considered to be a 32-bit real or 
integer value. When specific values are given in the table below, they are written to the binary file. 
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Number of words 


1 




ntents 


Version number of die input file = 5.0. 


Description of geometry read in record B 1 of section 2.2.1 
blank padded to 80 characters. Each character is written to 
MDDB.PLT as the real value of its integer ASCII value. 


Number of variables = 5.0 for which data is written to 
MDDBPLT. 


The 5-character names of the variables: 'NODE , 
'COLOR', 'X 'Y ’, 'Z ’. Each character is written 
to MDDB.PLT as the real value of its integer ASCII value. 


Ten character zone name: 'SURFACES . Each character 
is written to MDDB.PLT as the real value of its integer 
ASCII value. 


Format of data = 3.0, which indicates FEPOINT. 




Number of nodes NODETV. 


Number of elements in connectivity matrix NCONV. 


End of header marker = 


Zone marker = 299. 


Number of repeat variables = 0.0. 


One word of data for each of five variable values for the 
NODETV nodes. Data are written in order as 
NODE1, COLOR1, XI, Yl, Z1 
NODE2, COLOR2, X2, Y2, Z2 

etc. 


One word of data for each of the indices of the four points 
of the connectivity matrix for the NCONV elements. 
When the four nodes are connected in order, they form a 
quadrilateral (or a triangle if two adjacent points are 
identical). 


TECPLOT ASCH Output File from MDDB. MDDB also produces an ASCII 
TECPLOT file whenever a TECPLOT display of the object is selected. This file is readable by 
PREPLOT which creates a binary file used by TECPLOT. The ASCII file is named MDDB.PRE 
and is opened in subroutine GENSURF with a FORTRAN open statement as unit 7, STATUS = 
UNKNOWN, FORM = TORMATTED’. STATUS = UNKNOWN causes any existing 
version of MDDB.PRE to be overwritten. The new version of the file is written in subroutine 
CONNECM. The format of this file is described below. When specific values are given, they are 
written to the file. 


Record FI Format (1X, I TITLE=”’,A,"") 


Entry 1 Columns 

Type _ 


Description of geometry from record B 1 | 9-69 

character 

HEADER 1 


Record F2 Format (' VARIABLES=NODE,COLOR,X,Y,Z') 


B3 

























This record provides variable names for the columns of data in record F4. 


Record F3 Format (' ZONE T="SURF 1' 

\I=\I4,'J=\I4,T=FEP< 

DINT') 



Type 


Number of nodes (points) on the surfaces of the 
structure plus the surface normal direction arrows. 

20-23 

integer 

NODETV 

Number of elements in the connectivity matrix. 

WE5BS3H ! i 




Record F4 Format (I6,F10.4,3F10.3)I6,F10.4,3F10.3)') L,COLOR,PT 

This record is repeated NODETV times to generate all of die surfaces and surface normal 
direction arrows. 



Columns 

1 T ype 

Variable 

Node number. 

" mr ' 1-6 


L 

Color of node. 

7-16 

real 

■■gszsm 

X coordinate of node 

17-26 

real 

PT(1) or 
ARROW(U) 

Y coordinate of node 

mm 

real 


Z coordinate of node 

37-46 

real 

WM33SEMM 

WmBSfi ¥m 


Record F5 Format (4110) 

This record is repeated NCONV times, once for each element of the connectivity matrix. 



Columns 

—Type 

■kb 

The indices of the four nodes defining an element 
of the connectivity matrix. When the four nodes 
are connected in order, they form a quadrilateral (or 
a triangle if two adjacent points are identical). 

1-10, 11-20, 
21-30, 31-40 

integer 

m 
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